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1 

Menetelma vastaanottimen tahdistamiseksi, paikannusjarjestelma, 
vastaanotin ja elektroniikkalaite 

Nyt esilla oleva keksinto kohdistuu oheisen patenttivaatimuksen 1 joh- 
5 danto-osan mukaiseen menetelmaan vastaanottimen tahdistustoiminto- 
jen suorittamiseksi, oheisen patenttivaatimuksen 9 johdanto-osan mu- 
kaiseen paikannusjarjestelmaan, oheisen patenttivaatimuksen 14 joh- 
danto-osan mukaiseen vastaanottimeen, oheisen patenttivaatimuksen 
v 22 johdanto-osan mukaiseen elektroniikkalaltteeseen, seka oheisen 

10 patenttivaatimuksen 29 johdanto-osan mukaiseen elektroniikkalaittee- 
seen. 

Eras tunnettu paikannusjarjestelma on GPS-jarjestelma (Global Po- 
sitioning System), joka kasittaa talla hetkella yli 20 satelliittia, joista sa- 

15 manaikaisesti vastaanottimen nakyvissa on 4 tai enemman. esim. 
Suomessa voi havaita leveyspiirista riippuen jopa yli 14 satelliittia sa- 
manaikaisesti johtuen nakyvyydesta yli pohjoisnavan. Nama satelliitit 
lahettavat mm. satelliitin ratatietoa (Ephemeris data) seka tietoa satellii- 
tin kellonajasta. Sijainnin maarityksessa kaytettava vastaanotin paatte- 

20 lee sijaintlnsa normaalisti slten, etta vastaanottimessa lasketaan use- 
ammasta paikannusjarjestelmaan kuuluvasta satelliitista samanaikai- 
sesti lahetettavan signaalin kulkuaika vastaanottimeen. Sijainnin maa- 
rittamiseksi on vastaanottimen vastaanotettava tyypillisesti vahintaan 
neljan nakyvissa olevan satelliitin signaali, jotta sijainti voidaan laskea. 

25 

Jokainen GPS-jarjestelm§n toimiva satelliitti lahettaa ns. L1-signaalia 
1575,42 Mhz:n kantoaaltotaajuudeila. Tata taajuutta merkitaan myos 
154f 0 , missa f 0 -10,23 MHz. Lisaksi satelliitit lahettavat L2-signaalia 
1227,6 Mhz:n kantoaaltotaajuudeila, eli I20f 0 . Satelliitissa suoritetaan 

30 naiden signaalien modulointi ainakin yhdelia valesatunnaissekvenssllla. 
Kullakin satelliitilla tama valesatunnaissekvenssi on erilainen. Modu- 
loinnin tuloksena muodostuu koodimoduloitu laajakaistasignaali. Kay- 
tetty modulointitekniikka mahdollistaa sen, etta vastaanottimessa pysty- 
taan erottamaan eri satelliittien lahettamat signaalit, vaikka lahetyk- 

35 sessa kaytettavat kantoaaltotaajuudet ovat olennaisesti samat. Tasta 
modulointitekniikasta kaytetaan nimitysta koodijako-monikayttotekniikka 
(CDMA, Code Division Multiple Access). Kussakin satelliitissa L1 -sig- 
naalin moduloinnissa kaytetaan valesatunnaissekvenssina mm. ns. 
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C/A-koodia (Coarse/Acquisition code), jona kaytetaan Gold-koodia. Jo- 
kainen GPS-satelliitti lahettaa signaalia kayttamalla yksilollista C/A- 
koodia. Koodit muodostetaan kahden 1 023-bittisen binSarisekvenssin 
modulo-2 summana. Ensimmainen binaarisekvenssi G1 on muodos- 
5 tettu poiynomilla X 10 +X 3 +1 ja toinen binaarisekvenssi G2 on muodos- 
tettu viivastamalla polynomia X 10 +X 9 +X a +X 6 +X 3 +X 2 +1 siten, etta kulla- 
kin satelliitilla viive on erilainen. Tama jarjestely mahdollistaa sen, etta 
eri C/A-koodit voidaan muodostaa samanlaisella koodigeneraattorilla. 
C/A-koodit ovat siis binaarikoodeja, joiden kellotusnopeus (Chipping 
10 rate) GPS-jarjestelmassa on 1,023 MHz. C/A-koodi kasittaa 1023 alib'rt- 
tia (Chip), jolloin koodin toistoaika (koodijakso, epoch) on 1 ms. L1-sig- 
naalin kantoaaltoa moduloidaan viela navigointi-informaatiolla 50 bit/s 
bittinopeudella. Navigointi-informaatio kasittaa tietoa satelliitin 
"terveydentilasta" (health), radasta, aikatietoa jne. 

15 

Satelliitit tarkkailevat laitteistonsa kuntoa toimintansa aikana. Sateilittit 
voivat kayttaa esim, ns. vahtikoiratoimintoja joidenkin laitteistoon mah- 
dollisesti tulleiden vikojen havaitsemiseen ja ilmoittamiseen. Virheet ja 
toimintahairiot voivat oila hetkellisia tai pidempiaikaisia. Terveydentila- 

20 tietojen perusteella voidaan mahdollisesti osa virheista kompensoida tai 
jattaa vikaantuneen satelliitin lahettama informaatio kokonaan huomioi- 
matta. Lisaksi tilanteessa, jossa useamman kuin neljan satelliitin sig- 
naali on vastaanotettavissa, voidaan terveydentilatietojen perusteella 
painottaa eri satelliiteista vastaanotettua informaatiota eri tavalla, Tal- 

25 loin epaluotettavalta vaikuttavien satelliittien mahdollisesti aiheuttamia 
virheita mittauksiin voidaan minimoida. 

Satelliittien signaalien havaitsemiseksi ja satelliittien tunnistamiseksi on 
vastaanottimen suoritettava tahdistus, jossa vastaanotin etsii kulloinkin 
30 kunkin satelliitin signaalin ja pyrkii tahdistumaan tahan signaaliin, jotta 
signaalin mukana lahetettava data voidaan vastaanottaa ja demodu- 
loida. 

Sijainninmaaritysvastaanottimen on suoritettava tahdistus mm. silloin, 
35 kun vastaanotin kytketaan paalle ja myos tilanteessa, jossa vastaanotin 
ei ole pitkaan aikaan pystynyt vastaanottamaan minkaan satelliitin sig- 
naalia. Mm. kannettavissa laitteissa tallainen tilanne voi syntya helposti, 
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koska laite liikkuu ja laitteen antenni ei aina ole opfimaalisessa asen- 
nossa satelliitteihin nahden, mika heikentaa vastaanottimeen tulevan 
signaalin voimakkuutta. Myos kaupunkialueilla rakennukset, esim. nk. 
blokkaus eli signaalitien katkaisu, vaikuttavat vastaanotettavaan sig- 
5 naaliin ja lisaksi voi syntya ns. monitie-etenemista, jossa lahetetty sig- 
naali saapuu vastaanottimeen eri kulkureitteja, esim. suoraan satelliitis- 
ta (line-of-sight) ja lisaksi rakennuksista heijastuneena. Tama monitie- 
eteneminen aiheuttaa sen, etta sama signaali vastaanotetaan useina 
eri vaiheisina signaaleina. 

10 

Sijainninmaaritysjarjestelylle on kaksi paaasiallista tehtavaa: 

1. vastaanottimen pseudo-etaisyyden laskenta eri GPS-satel- 

liitteihin, ja 

15 2. vastaanottimen sijainnin maaritys, jossa kaytetaan laskettuja 

pseudo-etaisyyksia seka satelliittien sijaintitietoa. Sateliiittien 
kulloinenkin sijaintitieto voidaan laskea satelliiteista vas- 
taanotettujen ephemeris- ja aikakorjaustietojen perusteella. 

20 Etaisyyksia satelliitteihin nimitetaan pseudo-etaisyyksiksi, koska aika ei 
vastaanottimessa ole tarkasti tiedossa. Koska aikaa ei tiedeta abso- 
luuttisen tarkasti, on paikka ja aika selvitettava sopivimmin iteroimalla 
mitattua dataa linearisoidulla yhtaloryhmalla. Talloin sijainnin ja ajan 
maaritysta toistetaan, kunnes on saavutettu riittava tarkkuus ajan ja si- 

25 jainnin suhfeen. 

Pseudo-etaisyyden laskenta voidaan suorittaa mittaamalla eri satelliit- 
tien signaalien keskinaiset, naennaiset kulkuviiveet. Sen jalkeen, kun 
vastaanotin on tahdistunut vastaanotettuun signaaliin, suoritetaan sig- 
30 naalissa lahetetyn informaation llmaisu. 

Lanes kaikki tunnetut GPS-sijainninrnaaritysvastaanottimet (GPS-navi- 
gointivastaanottimet) kayttavat korrelaatiomenetelmia tahdistustoiminto- 
jen suorittamiseen, eli koodiin tahdistumiseen (acquisition) ja/tai tahdis- 
35 tuksen yllapitamiseen (seurantaan, tracking). Sijainninmaaritysvas- 
taanottimessa on tallennettu tai generoidaan paikallisesti vertailukoodit 
r(x), eli eri satelliittien valesatunnaissekvenssit. Vastaanotetulle sig- 
naalille suoritetaan muunto valitaajuudelle (Down Conversion), minka 
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jalkeen vastaanotin suorittaa vastaanotetun signaalin kertomisen tal- 
lennetulla valesatunnaissekvenssilla. Kertolaskun tuloksena muodostu- 
nut signaali integroidaan tai alipaastosuodatetaan, jolloin tuloksena 
saadaan tieto siita, onko vastaanotetussa signaalissa ollut jonkin satel- 
5 liitin lahettama signaali. Vastaanottimessa suoritettava kertolasku tois- 
tetaan siten, etta kullakin kerralla siirretaan vastaanottimeen tallennetun 
valesatunnaissekvenssin vaihetta. Oikea vaihe paatellaan korrelaatiotu- 
loksesta edullisesti siten, etta korrelaatiotuloksen ollessa suurin f on oi- 
kea vaihe loytynyt. Talloin vastaanotin on oikein tahdistunut vastaan- 
10 otettuun signaaliin. 

Sen jalkeen, kun koodiin tahdistuminen on suoritettu, suoritetaan viela 
taajuuden hienosaato ja vaihelukitus. Tama korrelaatiotulos ilmaisee 
myos GPS-signaalissa lahetetyn informaation. 

15 

Edella mainittu tahdistus ja taajuudensaatoprosessi on toistettava kul- 
lekin sellaisen satelliitin signaalille, jota vastaanottimessa vastaanote- 
taan. Tahan prosessiin kuluu siis runsaasti aikaa erityisesti tilanteessa, 
jossa vastaanotettavat signaalit ovat heikkoja. Taman prosessin no- 

20 peuttamiseksi kaytetaan joissakin tunnetun tekniikan mukaisissa vas- 
taanottimissa useampia korrelaattoreita, jolloin samanaikaisesti voidaan 
etsia useampia korrelaatiohuippuja. Pelkastaan korrelaattoreiden luku- 
maaraa kasvattamalla ei kaytannon sovelluksissa tahdistus- ja taajuu- 
densaatdprosessia voida nopeuttaa kovinkaan paljon f koska korrelaat- 

25 toreiden lukumaaraa ei vol kasvattaa rajattomasti. 

Joissakin tunnetun tekniikan mukaisissa GPS-vastaanottimissa on kay- 
tetty FFT-tekniikkaa vastaanotetun GPS-signaalin Doppler-siirtyman 
maarittamiseksi tavallisten korrelaattoreiden yhteydessa. Naissa vas- 
30 taanottimissa kaytetaan korrelointia vastaanotetun signaalin kaistanle- 
veyden kaventamiseksi 10 kHz-30 kHziiin. Tama kapeakaistainen sig- 
naali analysoidaan FFT-algoritmeilla kantoaaltotaajuuden maarittami- 
seksi. 

35 Kansainvalisessa patenttihakemuksessa WO 97/14057 on esitetty 
GPS-vastaanotin ja menetelma GPS-signaalien kasittelemiseksi. Tassa 
julkaisussa esitetty vastaanotin kasittaa paaasiassa kaksi erillista vas- 
taanotinta, joista ensimmainen vastaanotin on tarkoitettu kaytettavaksi 
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tilanteessa, jossa vastaanotetun signaalin voimakkuus on riittavan suuri 
ja toista vastaanotinta on tarkoitettu kaytettavaksi tilanteessa, jossa 
vastaanotetun signaalin voimakkuus ei riita riittavan tarkkaan sijainnin 
maaritykseen kaytettaessa ensimmaista vastaanotinta. Tassa toisessa 
5 vastaanottimessa suoritetaan vastaanotetun signaalin digitointi ja tal- 
lennus muistivalineisiin, jolloin myohemmin nSita tallennettuja signaalei- 
ta kasitellaan digitaalisessa signaalinkasittely-yksikossa. Digitaalinen 
signaalinkasittely-yksikko suorittaa konvoiuutio-operaatioita vastaanote- 
tulle digitoidulle signaalille. Naiden konvoluutio-operaatioiden tavoittee- 
10 na on pseudo-etaisyyksien laskenta. Muistivalineisiin tallennetaan tyy- 
pillisesti sadasta tuhanteen koodijaksoa (PM Frames), mika vastaa 
100 ms-1 s mittaista signaalia. Taman jalkeen tutkittavan satelliitin 
koodia vastaava tallennettu koodi noudetaan vastaanottimen muistista 
kaytettavaksi vastaanotetun signaalin analysointiin. 

15 

Vastaanottimessa poistetaan myos Doppler-siirtyma. Taman Doppler- 
siirtyman maara on selvitetty joko ensimmaisella vastaanottimella tai 
GPS-jarjestelmaan kuuluvasta tukiasemasta vastaanotetun tiedon pe- 
rusteella. Taman jalkeen suoritetaan perakkaisten kehysten koherentti 

20 summaus, Talle summauksen tuloksena saadulie datajoukoile suorite- 
taan nopea Fourier-muunnos. Fourier-muunnostulokselle suoritetaan 
kertolasku, jossa toisena kertojana kaytetaan muistivalineisiin tallenne- 
tun referenssisignaalin Fourier-muunnoksen kompleksikonjugaattia, 
Talle kertolaskun tulokselle suoritetaan edelleen kaanteinen Fourier- 

25 muunnos, jolloin saadaan joukko korrelaatiotuloksia. Tassa julkaisussa 
on siis korrelaatio korvattu Fourier-muunnoksella, jolloin laskutoimitus- 
ten maaraa on saatu pienennettya. Julkaisun mukaan menetelma no- 
peuttaa 10-100 -kertaisesti sijainninmaSritysta verrattuna mainitun jul- 
kaisun hakemishetkella tunnettuihin ratkaisuihin. 

30 

Joissakin tunnetun tekniikan mukaisissa vastaanottimissa kaytetaan 
sovitettuja suodattimia (matched filter) signaaliin lukittumisessa ja seu- 
rannassa. Tallaisten sovitettujen suodattimien kaytossa on kuitenkin se 
epakohta, etta ne kuluttavat suhteellisen paljon tehoa, mika rajoittaa 
35 niiden kayttoa erityisesti kannettavissa laitteissa. 
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Nyt esilla olevan keksinnon eraana tarkoituksena on aikaansaada vas- 
taanotin, jossa tahdistus lahetettyyn signaaliin voidaan suorittaa merkit- 
tavassa maarin nopeammin kuin tunnetun tekniikan mukaisissa vas- 

5 taanottimissa. Keksinto soveltuu erityisesti sijainninmaaritysvastaan- 
ottimissa kaytettavaksi, mutta rnyos mulssa vastaanottimissa, edulli- 
sesti CDMA-vastaanottimissa, joissa vastaanottimen on tahdistuttava 
hajaspektrisignaaliin. Keksinto perustuu siihen ajatukseen, etta suorite- 
taan aika/taajuus-muunnokset seka kaanteis-, eli taa- 

10 juus/aikamuunnokset nopealla Hartley-muunnoksella ja vastaavalla 
kaanteismuunnoksella. Talloin voidaan vastaavat laskutoimitukset 
suorittaa nopeammin kuin tunnetun tekniikan mukaisissa ratkaisuissa, 
kuten FFT-muunnoksiin perustuvissa jfirjestelmissa. Nyt esilla olevan 
keksinnon mukaiselle menetelmalle on tunnusomaista se, mita on esi- 

15 tetty oheisen patenttivaatimuksen 1 tunnusmerkkiosassa. Nyt esilla 
olevan keksinnon mukaiselle paikannusjarjestelmalle on tunnusomaista 
se, mita on esitetty oheisen patenttivaatimuksen 9 tunnusmerkkiosas- 
sa. Nyt esilla olevan keksinnon ensimmalsen edullisen suoritusmuodon 
mukaiselle vastaanottimelle on tunnusomaista se f mita on esitetty ohei- 

20 sen patenttivaatimuksen 14 tunnusmerkkiosassa, Nyt esilla olevan 
keksinnon mukaiselle elektroniikkalaitteelle on tunnusomaista se, mita 
on esitetty oheisen patenttivaatimuksen 22 tunnusmerkkiosassa. Nyt 
esilla olevan keksinnon toisen edullisen suoritusmuodon mukaiselle 
elektroniikkalaitteelle on viela tunnusomaista se, mita on esitetty ohei- 

25 sen patenttivaatimuksen 29 tunnusmerkkiosassa. 

Nyt esilla olevalla keksinnolla saavutetaan merkittavia etuja tunnetun 
tekniikan mukaisiin menetelmiin ja vastaanottimiin verrattuna. Keksin- 
non mukaisella menetelmalla voidaan muunnokset suorittaa nopeam- 

30 , mln ja pienemmalla tehonkulutuksella verrattuna esim. FFT-pohjaisiin 
jarjestelmiin. Talloin erityisesti sijainninmaaritysvastaanottimissa voi- 
daan vastaanottimen tahdistusta vastaanotettavaan signaaliin nopeut- 
taa. Lisaksi keksinnon mukaisen menetelman mukainen reaaliaikainen 
aika/taajuusmuunnos ja taajuus/aikamuunnos toteuttaa pinta-alaltaan 

35 merkittavasti pienemmalla puolijohdepiirilla kuin tunnetun tekniikan mu- 
kaiset ratkaisut. Keksinnon mukaisella menetelmalla voidaan saavuttaa 
merkittava parannus tahdistusnopeuteen verrattuna tunnetun tekniikan 
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mukaisiin menetelmiin. Keksinnon mukaisen vastaanottimen kokonais- 
energiankulutus voidaan pitaa kohtuullisena, jolloin keksinto soveltuu 
erityisen hyvin kaytettavaksi kannettavissa laitteissa. Talloin sijainnin- 
maaritysvastaanotin voidaan toteuttaa myos langattoman viestimen 
5 yhteyteen. 

Nyt esilla olevaa keksintoa selostetaan seuraavassa tarkemmin viitaten 
samalla oheisiin piirustuksiin, joissa 



10 kuva 1 esittaa pelkistettyna. lohkokaaviona vastaanotinta, jossa 
keksinnon mukaista menetelmaa voidaan soveltaa, 



15 



kuva 2 esittaa keksinnOn era&n edullisen suoritusmuodon mukaisen 
menetelman korreloinfivaihetta pelkistettyna kaaviona, 

kuva 3 esittaa keksinnon eraan edullisen suoritusmuodon mukaisen 
menetelman analyysivaihetta pelkistettyna kaaviona, 



kuva 4 esittaa keksinnSn eraan edullisen suoritusmuodon mukaisen 
20 menetelman paattelyvaihetta pelkistetysti, 

kuva 5 esittaa keksinnon eraan edullisen suoritusmuodon mukaista 
elektroniikkalaitetta pelkistettya lohkokaaviona, ja 

25 kuva 6 esittaa yksinkertaistettuna periaatekaaviona keksinnon 
eraan edullisen suoritusmuodon mukaista paikannusjarjes- 
telmaa. 



Seuraavassa keksintoa selostetaan kayttaen esimerkkina kuvan 1 mu- 
30 kaista GPS-vastaanotinta. Tarkastellaan aluksi menetelmassa sovellet- 
tavien Hartley-muunnoksen ja kaanteisen Hartley-muunnoksen ominai- 
suuksia. Hartley-muunnos voidaan esittaa kaavalla: 

KO = S *(^o ) • [cos(27tfkt o / M) + sin{2nfkt 0 /M)]/M (1 ) 

35 

ja kaanteinen Hartley-muunnos voidaan esittaa kaavalla: 
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= E H { M o ) ■ [cos(2tc/c/ 0 1 1 M) + s\n(2nkf 0 1 1 M)] (2) 
missa 

5 x(k) ovat mitatut naytteet aikatasossa, 

M on naytteiden lukumaara, 
to on naytteenottovali, ja 
fo on 1/[(M-1)to] 

10 Diskreetin Hartley-muunnoksen laskemiseksi tarvitaan edella esltettyjen 
kaavojen mukaan N 2 laskutoimitusta. Tata maaraa voidaan kuitenkin 
pienentaa matemaattisten muunnosten teoriasta tunnettujen permutaa- 
tioiden ja ns. Butterfly-prosessin avulla, jolloin muunnokset voidaan 
suortttaa Nlog 2 N laskutoimituksella. Tata muunnosta kutsutaan talldin 

15 nopeaksi Hartley-muunnokseksi FHT (Fast Hartley-transform). Hartley- 
muunnoksen ominaisuuksien perusteella konvoluutio aikatasossa voi- 
daan suorittaa kerto- ja yhteenlaskujen avulla taajuustasossa seuraa- 
vasti: 

20 x,(l)** 2 (f)-H,(0- ^aW + H|(-0-«feoW (3) 
missa 

Xi(t) ja x 2 (t) ovat konvoloitavat signaalit, 
Hi(f) ja H 2 (f) ovat x,(t) ja X2(t) Hartley-muunnokset vastaa- 
25 vasti, ja 

H 2e (f) ja H 2o (f) ovat signaalln X2(t) Hartley-muunnoksen paril- 
linen ja pariton muunnos vastaavasti. 

Jos ainakin yksi konvoloitava signaali on pariton tai parillinen, seuraa 
30 siita se, etta konvoluutioyhtalo yksinkertaistuu muotoon H,(f) H 2 (f). 

Edella esitettya konvoluutiota voidaan kayttaa hyodyksi mm. lasket- 
taessa kahden signaalin valista korrelaatiota. Tailoin korrelaation las- 
kemiseksi konvoluution avulla, kaannetaan jompi kumpi signaaleista ai- 
katasossa peilikuvaksi ja lasketaan konvoluutio edullisesti edella esitet- 
35 tyja Hartley-mLiunnoksia kayttamalla. Tata on kuvattu jaljempana tassa 
selityksessa. 
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Kuvataan seuraavaksi keksinnon ensimmaisen edullisen suoritusmuo- 
don mukaisen menetelman vastaanottimen 1 toimintaa viitaten samalla 
kuvaan 1. Vastaanottimessa 1 vastaanotettava signaali muunnetaan 
5 sopivimmin valitaajuudelle muunninlohkossa 2. Tassa vaiheessa sig- 
naali kasittaa sinansa tunnetusti kaksi komponenttia: I- ja Q-kompo- 
nentit, joiden valilla on 90° vaihe-ero. Nama valitaajuudelle muunnetut 
anaiogiset signaalikomponentit digitoldaan digitointilohkossa 3 ja johde- 
taan kertojalohkoon 4. Kertojalohkossa 4 digitoidun signaalin I- ja Q- 

10 komponentit kerrotaan numeerisesti ohjatun oskillaattorin 5 (NCO, Nu- 
merically Controlled Oscillator) muodostamalla signaalilla. Tama nu- 
meerisesti ohjatun oskillaattorin 5 signaali on tarkoitettu korjaamaan 
Doppler-siirtymasta ja vastaanottimen 1 paikallisoskillaattorin (el esltet- 
ty) taajuusvirheesta johtuva taajuuspoikkeama. Kertojalohkon 4 muo- 

15 dostama signaali johdetaan tahdistuslohkoon 6, jossa keksinnon mu- 
kaista menetelmaa sovelletaan. Tama tahdistuslohko 6 pyrkii siis loy- 
tamaan satelliitin koodivaiheen ja taajuuspoikkeaman kaytettavaksi 
tahdistuksen jalkeisissa toiminnoissa. Tata kuvataan jaljempana tassa 
selityksessa. Ohjauslohkolla 7 ohjataan tahdistuksen aikana skannaus- 

20 lohkoa 8, jonka avulla numeerisesti ohjatun oskillaattorin 5 taajuutta 
saadetaan tarvittaessa. Ohjauslohko 7 ohjaa ensimmaista kytkinta 9 
kytkemaan joko taman skannauslohkon 8 muodostaman signaalin tah- 
distuksen aikana numeerisesti ohjatulle oskillaattorille 5 tai koodin selvi- 
tyslohkon 1 1 muodostaman ohjaussignaalin numeerisesti ohjatulle os- 

25 killaattorille 5 sen jalkeen, kun tahdistus on saavutettu. Toisella kytki- 
mella 10 ohjataan tata seurantalohkon 11 toimintaa. Tama seuranta- 
lohko 11 muodostaa osan sinansa tunnetusta koodivaihelukitusta sil- 
mukasta ja kantoaaltovaihelukitusta silmukasta (ei esitetty). 

30 Kayttojannitteiden kytkemisen jalkeen tal tllanteessa, jossa vastaanotin 
1 ei ole pystynyt vastaanottamaan GPS-satelliittien signaalia pitkaan 
aikaan, suoritetaan vastaanottimessa 1 edullisesti kaksidimensioinen 
etsintavaihe kullekin satelliitille, joiden signaalia vastaanotetaan. Tassa 
kaksidimensioisessa etsintavaiheessa tarkoituksena on selvittaa kunkin 

35 satelliitin kantoaaltotaajuus seka koodivaihe. Tahan kantoaaltotaajuu- 
teen vaikuttaa siis satelliitin liikkumisesta johtuva Doppler-siirtyma, 
monitie-eteneminen seka vastaanottimen paikallisoskillaattorin epatark- 
kuudet. Taajuusepavarmuus voi olla varsin suuri, jopa ± 6 kHz, jolloin 
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vastaanottimen 1 on suoritejtava etsinta n. 12kHz:n taajuusalueelta 
varsinaisen lahetystaajuuden 1 (L1 = 1575,42 MHz) suhteen. Vastaanotin 
1 ei myoskaan tieda tarkkaa koodivaihetta, jolloin vastaanottimen on 
suorltettava myos koodivaiheen selvitys 1023 mahdollisesta koodivai- 
5 heesta. Talloin saadaan kaksidimensioinen etsintaprosessi, jossa yksi 
dimensio on taajuuspoikkeama 12 kHz:n alueella ja toinen dimensio on 
koodivaihe 1023:sta erilaisesta koodivaiheesta. Keksinnon eraan edul- 
lisen suoritusmuodon mukaisessa menetelmassa voidaan kerrallaan 
tutkia n, 500 Hz:n taajuusalue, jolloin menetelmaa toistetaan tarvittaes- 

10 sa 24 kertaa koko tutkittavan 12 kHz:n taajuusalueen kattamiseksl. On 
selvaa, etta tassa selityksessa kaytettavat esimerkkian/ot ovat vain 
keksintoa selventavina, mutta eivat rajoittavina esimerkkeina. Keksln- 
toa voidaan sovoltaa muissakin kuin GPS-jarjestelmissa, jolloin mm. 
mainitut taajuusarvot, koodivaiheet seka koodien lukumaara voivat 

15 vaihdella. 

Selostetaan seuraavaksi keksinnon ensimmmaisen edullisen suoritus- 
muodon mukaisen menetelman toimintaa kuvan 1 mukaisessa vas- 
taanottimessa 1. Tahdistuksen kaynnistamiseksi skannauslohko 8 

20 asettaa numeerisesti ohjatun: osklllaattorin 5 taajuuden siten, etta vas- 
taanottimella 1 vastaanotetaan edullisesti taajuusalueen pienimpia 
taajuuksia, tassa esimerklssa 1575,414 MHz — 1575,4145 MHz. Vas- 
taanotin voi maarittaa aloitustaajuuden myos siten, etta vastaanotin 
kayttaa hyvaksi esim. aikaisemmin selvitettya sijaintitietoa ja/tai alma- 

25 nakkatietoa, jolloin sijainninmaaritysta voidaan viela nopeuttaa. Vas- 
taanotettavasta signaaJista tallennetaan naiytteita naytevektoreiden 
muodostuslohkossa 12 edullisesti kompleksisiksi naytevektoreiksi 
p(1), p(2)...p(N), joissa kussakin on 1023 naytetta tassa edullisessa 
suoritusmuodossa. Naytevektoreiden muodostuslohkon 12 naytteiden 

30 tallennustaajuus tassa edullisessa suoritusmuodossa on olennaisesti 
sama kuin alibittien (Chips) kellotusnopeus, eli n. 1 023 000 naytetta 
sekunnissa. Naytevektorit ovat jatkuvia siten, etta seuraava naytevek- 
tori jatkuu ajallisesti edellisen naytevektorin jalkeen, eli naytevektorin 
viimeisen naytteen ja seuraavan naytevektorin ensimmaisen naytteen 

35 aikaero on olennaisesti sama kuin naytevektorin perakkaisten nayttei- 
den valinen aikaero. Nama 1023 naytetta vastaavat siis n. 1 ms:n mit- 
taista signaalia, jolloin aika-taajuusmuunnoksessa taajuusalue on 
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n. 1 kHz, josta osa voidaan hyddyntaa. Naytevektoreiden muodostus- 
vaihetta on esitetty viitteella 1 01 oheisessa kuvassa 2. 

Naytevektoreita on edullisesti N kappaletta. Lisaksi naytevektoreiden 
5 muodostusta voidaan toistaa. Naytevektoreiden p(1), p(2)...p(N) luku- 
maaran N arvoa maaritettaessa tulee huomioida GPS-jarjestelmassa 
se, etta signaaliin on moduloitu informaatiota bittinopeudella 50 bittia/s 
binaarisena vaihemodulaationa. Toinen tata naytevektoreiden 
p(1), p(2)...p(N) lukumaaraa N rajoittava tekija on vastaanottimen pai- 
10 kallisoskillaattorin taajuusstabiilius. 

Naytevektoreiden muodostusvaiheen lisaksi keksinnon mukaisessa 
tahdistusmenetelmassa suoritetaan korrelointivaihe korrelaatiofunktio- 
matriisin muodostamiseksi. 

15 

Tama korrelointivaihe voidaan suorittaa osittain jo naytteenoton aikana. 
tai sen jalkeen kun N kappaletta naytevektoreita p(1), p(2)...p(N) on 
muodostettu. Jos korrelointivaihe suoritetaan esim. siten, etta kunkin 
naytevektorin tallennuksen jalkeen lasketaan sille aika-taajuusmuunnos 

20 Hartley-muunnoksella, sopivimmin nopealla Hartley-muunnoksella 
(FHT), voidaan samaa aika-taajuusmuunninta kayttaa kaikissa N kap- 
paleessa naytevektoreita p(1), p(2)...p(N). Jos sen sijaan korrelointivai- 
he suoritetaan N naytevektorin tallennuksen jalkeen, on kaytettava joko 
kullekin naytevektorille omaa aika-taajuusmuunninta, tai aika-taajuus- 

25 muunnokset suoritetaan eri naytevektoreille perakkain samassa aika- 
taajuusmuuntimessa. Kuvassa 2 on esitetty menetelman korrelointivai- 
hetta, jossa naytevektoreista p(1), p(2)...p(N) muodostetaan korrelaatio- 
funktiomatriisi C x . 

30 Kullekin naytevektorille p(1), p(2)...p(N) suoritetaan diskreetti Hartley- 
muunnos 102, sopivimmin nopea Hartley-muunnos, eli FHT-muunnos: 

P(i) = FHT(p(/}), missa 1=1,..., N (4) 

35 Tata esittavat lohkot FHT1, FHT2, .... FHTN kuvassa 2. Kaytanndn las- 
kutoimituksissa kaytetaan edullisesti 1024:aa arvoa, koska talloin 
diskreetti Hartley -muunnos on kaytanndn sovelluksissa huomattavasti 
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tehokkaammin toteutettavissa (FHT-algoritmilla) kuin kaytettaessa 
1023:a arvoa. Eras tapa on lisata ylimaarainen nolla 1024:ksi alkioksi. 
Talla on vahainen vaikutus muunnostulokseen. 

5 Vastaanottimessa on tallennettuna sopivimmin kaikkien GPS-jarjestel- 
man satelliittien (ei esitetty) C/A-koodia vastaavat vertailukoodit r(x), 
missa x viittaa satelliitin tunnukseen, esim. x on valilla 1 — 30. Vertailu- 
koodeja ei valttamatta tarvitse tallentaa, vaan ne voidaan myos gene- 
roida vastaanottimessa. Korrelaatiovaiheessa valitaan tai generoidaan 
10 kulloinkin sen satelliitin vertailukoodi, jonka lahettamaan signaaliin 
vastaanotinta tahdistetaan. Vertailukoodi on kaannetty ajallisesti taka- 
perin. Talle takaperoiselle vertailukoodille, jota kuvassa 2 on merkitty 
viitteella r(x) , suoritetaan diskreetti Hartley-muunnos 103, sopivimmin 

nopea Hartley-muunnos: 
15 _ 

R( x)£=FHT (r(x)) (5) 

Takaperoinen vertailukoodi r(x) ja/tai sen FHT-muunnos on usein tal- 

lennettu jo etukateen vastaanottimen muistivalineisiin, mutta se voidaan 
20 muodostaa vertailukoodista r(x) tahdistuksen yhteydessa. 

Seuraavaksi korrelaatiovaiheessa suoritetaan kunkin naytevektorin p(i) 
Hartley-muunnostuloksen P(i) ja takaperoisen vertailukoodin r(x) Hart- 

ley-muunnoksen W(x) valinen kertolasku 104: 
25 _ 

M x (i) = R(x)*P(i) (6) 

Naille kertolaskujen tuloksille suoritetaan viela kaanteinen Hartley- 
muunnos 105 (IFHT), jolloin tuloksena saadaan vertailukoodin r(x) ja 
30 vastaanotetun signaalin ristikorrelaatio kaikilla mahdollisilla kokonaislu- 
kuviiveilla (1023 kpl). 



m x 0) = »FHTT(M x (/)) 



(7) 
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Tama perustuu siihen edella esitettyyn seikkaan, etta aikatason signaa- 
lien konvoluution Hartley-muunnos vastaa aikatason signaalien Hartley- 
muunnosten kertolaskua. Kun lisaksi kaytetaan takaperoista vertailu- 
koodia, voidaan Hartley-muunnoksella suorittaa nopea diskreettiaikai- 
5 nen korrelaatio. Korrelaatiofulos kasittaa tassa edullisessa esimerkissa 
talloin 1023 alkiota. Naista eri naytevektoreista p(i) muodostetuista risti- 
korrelaatiotuloksista m x (i) muodostetaan korrelaatiofunktiomatriisi C XI 

jossa rivien maara on naytevektoreiden lukumaara N. 

10 On selvaa, etta vertailukoodin ajallisen kaantamisen sijasta voidaan 
naytevektoreista p(i) muodostaa takaperoisia naytevektoreita p(/) , jol- 

loin edella esitetyissa laskutoimituksissa kaytetaan vertailukoodia r(x) 
suoraan ja takaperoisia naytevektoreita, ErSassa edullisessa suoritus- 
muodossa el tatvitse tehda kumpaakaan edella mainituista kaantami- 
15 sista, eli voidaan kayttaa vertailukoodia r(x) ja naytevektoreita p(i) sel- 
laisenaan. 

Korrelaatiovaiheessa muodostetun korrelaatiofunktiomatriisin C x rivit 
esittavat vastaanotetun signaalin ja vertailukoodin ristikorrelaatiota eri 
20 vaihe-eroilla yhden millisekunnin valein otettuna. Kaavana korrelaatio- 
funktiomatriisi voidaan talloin esittaa seuraavasti: 



m X(k (1) 
m Kik (2) 



(8) 



25 Seuraavassa vaiheessa eli analyysivaiheessa kaytetaan korrelaatio- 
funktiomatriisin C x transpoosia 106, jossa rivit esittavat aikatasossa 

signaalin naytteita siten kuin tunnetun tekniikan mukaisessa korrelaat- 
torissa. Kukin rivi vastaa tiettya koodivaihe-eroa vastaanotetun signaa- 
lin ja vertailukoodin valilla, Taman korrelaatiofunktiomatriisin C x 

30 transpoosin kullekin riville suoritetaan Hartley-muunnos 107 koherentin 
etsintamatriisin A* muodostamiseksi, joiloin voidaan suorittaa taajuus- 

analyysi todellisen taajuussiirtyman selvittamiseksi. Tata esittaa ohei- 
nen kuva 3. 
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A x =FHT(cJ) (9) 

Kaytannon sovelluksissa ei korrelaatiofunktiomatriisista tarvitse erik- 
seen muodostaa transponoitua matriisia, vaan tallennetun korrelaatio- 
5 funktiomatriisin C x alkiot luetaan muistista 16 (kuva 5) eri suunnassa, 
edullisesti sarakkeittain. 

Tahdistuslohkossa 6 voidaan tarvittaessa suorittaa viela epakoherentti 
summausvaihe signaalikohinasuhteen parantamiseksi. Epakoherentin 
10 summausvaiheen toteuttamiseksi toistetaan edella esitettyja naytevek- 
toreiden muodostusvaihetta, korrelointivaihetta ja analyysivaihetta K 
kertaa (ei esitetty). Tama toistojen lukumaara K valitaan edullisesti si- 
ten, etta signaali-kohinasuhdetta saadaan riittavasti parannettua, kui- 
tenkin kohtuullisessa ajassa. 

15 

Selvitysvaiheessa (acquisition) tutkitaan etsintamatriisin A x alkioiden 
a x (ij) arvoja ja pyritaan Idytamaan 110 sellainen arvo, joka ylittaa en- 

nalta asetetun kynnysarvon ja on selvasti muita arvoja suurempi. Jos 
tallainen arvo loytyy, iimaisee se koodivaihe-eron seka taajuuspoik- 

20 keaman, koska kyseessa on todennakoisesti satelliitin lahettama sig- 
naali. Mikali signaali ei ole satelliitin lahettamaa, vaan kohinaa tai 
muuta satunnaista hairiosignaalia, ei merkittavia korrelaatiohuippuja 
p'rtaisi muodostua. Koodivaihe-ero ilmenee taman suurimman arvon 
rivi-indeksista ja vastaavasti taajuuspoikkeama ilmenee sarake-indek- 

25 sista. Kuvan 4 esimerkissa on esitetty eras tilanne, jossa selvasti muita 
arvoja suurempi arvo on loytynyt. Jos sen sijaan etsintamatriisista A x ei 
loydy tallaista arvoa, elj tutkitulla taajuusalueella ei todennakoisesti 
vastaanotettu haettavan satelliitin lahettamaa signaalia, muutetaan 
tutkittavaa taajuuskaistaa ja suoritetaan edella esitetyt vaiheet etsinta- 

30 matriisin muodostamiseksi. Talla menetelmalla saadaan koko tutkittava 
6 kHz:n alue kaytya lapi toistamalla edella esitetyt vaiheet tarvittavan 
monta kertaa. 

Tarvittaessa voidaan edella esitetyt vaiheet toistaa koko tutkittavalle 
35 taajuuskaistalle ja tallentaa eri toistokerroilla muodostetut etsintamat- 
riisit tai vain mahdolliset huippukohdat, ennen kuin etsitaan suurinta 
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korrelaatiohuippua. Talloin voidaan pienentaa virhetulkintojen mahdolli- 
suutta mm. sellaisessa tilanteessa, jossa kynnysarvo on asetettu liian 
pieneksi ja hairiosignaali voi aiheuttaa virhetulkinnan, 

5 Sen jalkeen kun oikea taajuuspoikkeama ja koodivaihe on selvitetty, 
voidaan vastaanotin asettaa seurantatilaan (tracking). Heikoimmilla 
signaaleilla ei datavastaanotto onnistu, vaan on sinansa tunnetusti tur- 
vauduttava esim. matkaviestinverkon kautta saatuun dataan. Etai- 
syysmittaus on edelleen mahdollista alentuneella tarkkuudella, Seuran- 
10 tatila asetetaan vaihtamalla kytkimet (kuva 1) toiseen asentoon, jolloin 
vastaanotettu informaatio johdetaan seurantalohkoon 11, jossa muo- 
dostetaan myos takaisinkytkenta numeerisesti ohjatun oskillaattorin 5 
taajuuden hienosaatoa varten. 

15 Sijainnin laskemiseksi vastaanotin suorittaa signaalin vastaanoton so- 
pivimmin vahintaan neljasta satelliitista vastaanotetun signaalin perus- 
teella. Talloin edella esitetty tahdistus toistetaan tarvittaessa kunkin 
satelliitin signaalia varten, jolloin vertailusekvenssiksi r(x) valifaan kul- 
loinkin sen satelliitin koodi, johon tahdistus suoritetaan. 

20 

Suuri osa keksinnon ensimmaisen edullisen suoritusmuodon mukaisen 
menetelman toteuttamisessa tarvittavista lohkoista voidaan toteuttaa 
esim. digitaalisessa signaalinkasittely-yksikossa (Digital Signal Proces- 
sor, ei esitetty). FHT-muunnosten suorittamisessa voidaan kayttaa joko 
25 laitteistopohjaisia ratkaisuja tai digltaalisen signaalinkasittely-yksikon 
ohjelmallisia toteutuksia. Lisaksi vastaanottimen toiminnan ohjaamiseen 
voidaan kayttaa ohjauselinta, edullisesti mikroprosessoria tai vastaa- 
vaa. 

30 Kayttamalla Hartley-muunnosta aika-taajuus- ja taajuus-aikamuunnok- 
sena saavutetaan viela se etu, etta molemmissa muunnoksissa voi- 
daan kayttaa samaa algoritmia, mika voidaan havaita kaavoja (1) ja (2) 
vertailemalla. Laitteistopohjaisessa ratkaisussa muunnosten tarvitsema 
elektroniikka voidaan toteuttaa joko noin puolta pienemmalla piiripinta- 

35 alalia (esim. ASIOpiirilla) ja tehonkulutuksella kuin vastaava toteutus 
Fourier-muunnoksia kaytettaessS prosessointiajan ollessa samaa suu- 
ruusluokkaa, tai samansuuruisella piiripinta-alalla prosessointinopeu- 
den jopa kaksinkertaistuessa. 
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Hartley-muunnosten muunnostuloksia H(f) voidaan kayttaa viela sig- 
naalin amplitudispektrin ja tehospektrin laskemiseen. Amplitudispektri 
voidaan laskea kaavoilla: 

5 

x real = H{f) + f ) (reaaliosa) (1 0) 

Ximag - MO - H(N- f) (imaginaariosa) (1 1 ) 

Tehospektri on laskettavissa kaavalla: 

10 

S = [H(0 2 +H(A/-0 2 ]/2 (12) 

Oheisessa kuvassa 5 on esitetty viela keksinnon eraan edullisen suori- 
tusmuodon mukainen elektroniikkalaite 24, joka kasittaa langattoman 

15 viestimen ja sijainninmaaritysvastaanottimen. Ensimrnaista antennia 13 
kaytetaan paikannussatelliittien lahettaman signaalin vastaanottami- 
seen. Vastaanotettu signaali johdetaan ensimmaiseen radio-osaan 14, 
jossa signaali muunnetaan valitaajuudelle ja digitoidaan. Ensimmainen 
radio-osa kasittaa mm. kuvan 1 vastaanottimen muunninlohkon 2, digi- 

20 tointilohkon 3 ja kertojalohkon 4. Digitoitu signaali, joka tassa vaiheessa 
kasittaa edullisesti I- ja Q-komponentit, johdetaan digitaaliseen signaa- 
likasittely-yksikkoon 15, jossa suoritetaan mm. naytevektoreiden muo- 
dostus naytevektoreiden muodostuslohkossa 12. Naytteita tallennetaan 
ensimmaisiin muistivalineisiin 16, jotka kasittavat mm. 

25 luku/kirjortusmuistia seka edullisesti myos lukumuistia ja/tai haihtuma- 
tonta luku/kinoitusmuistia digitaalisen signaalinkasittely-yksikon 15 oh- 
jelmakoodin tallennusta varten. Tassa suoritusmuodossa on signaalin- 
kasittely-yksikossa 15 toteutettu mm. tahdistuslohkon 6 toimintoja, ku- 
ten korrelaatiofunktiomatriisin C xk muodostus esim. Hartley-muuntimilla 

30 FHT1, FHT2 FHTN ja/tai sovitetuilla suodattimilla. Digitaalisessa 

slgnaalinkasittely-yksikossa 15 muodostetaan edullisesti myos kohe- 
rentti etsintamatriisi A x ja suoritetaan etsintavaihe. Digitaalinen signaa- 
linkasittely-yksikko 15 valittaa tiedon lasketusta vaihe-erosta ja taajuus- 
poikkeamasta suoritinlohkoon 17, joka kasittaa esim. mikroprosessorin 

35 ja l/O-logiikkaa. Suoritinlohko 17 ohjaa skannauslohkoa 8 seka ensim- 
rnaista kytkinta 9. Seurantalohko 1 1 on edullisesti ainakin osittain to- 
teutettu suoritinlohkon ohjelmakaskyina. Suoritinlohkon 17 datamuistina 
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ja ohjelmamuistina kaytetaan toisia muistivalineita 18. On selvaa, etta 
ensimmaiset rnuistivalineet 16 ja toiset rnuistivalineet 18 voivat kasittaa 
myos yhteista muistia. Sijainninmaaritysinformaatlota voidaan esittaa 
kayttajalle naytolla 19, 

5 

Suoritinlohkon 17 sovellusohjelmistossa on toteutettu myos langatto- 
man viestimen toimintoja. Nayttoa 19 voidaan talloin kayttaa sinansa 
tunnetusti myos mm. puheluinformaation esittamiseen. NSppalmistolla 
20 kayttaja voi ohjata sijainnimaaritysvastaanotinta seka langatonta 
10 viestinta. Koodekin 21 avuila suoritetaan audiosignaalien koodaarninen 
ja dekoodaaminen. Kuvassa 5 on esitetty viela langattoman viestimen 
radio-osa 22 ja toinen antenni 23. 

Edella kuvatussa keksinnon ensimmaisen edullisen suoritusmuodon 
15 mukaisessa menetelmassa sijainninmaaritys suoritettiin vastaanottimen 
yhteydessa. Keksinnon eraan toisen edullisen suoritusmuodon mukai- 
sessa menetelmassa ainakin osa sijainninmaarityslaskennasta suorite- 
taan tiedonsiirtoverkon yhteydessa. Tiedonsiirtoverkkoon voi kuulua 
esimerkiksi matkaviestinverkko. Talloin tiedonsiirtoverkkoon tiedonsiir- 
20 toyhteyteen on jarjestetty laskentapalvelin 26 (kuva 6) tai vastaava si- 
jainninmaarityslaskentaan soveltuva yksikko. Eras mahdollinen lasken- 
tapalvelimen sijoituspaikka on matkaviestinverkon tukiasema 25, mutta 
myos muita vaihtoehtoisia sijoituspaikkoja voidaan soveltaa nyt esitla 
olevan keksinnon yhteydessa. 

25 

Selostetaan seuraavaksi keksinnon eraan toisen edullisen suoritus- 
muodon mukaisen menetelman toimintaa kuvan 5 mukaisessa elek- 
troniikkalartteessa ja kuvan 6 mukaisessa jarjestelmassa. Tassa mene- 
telmassa suoritetaan olennaisesti samoja toimenpiteita kuin keksinnon 

30 ensimmaisen edullisen suoritusmuodon mukaisessa vastaanottimessa 
1, joten seuraavassa on kuvattu vain menetelman paapiirteita. Tahdis- 
tuksen kaynnistamiseksi skannauslohko 8 asettaa numeerisesti ohjatun 
oskillaattorin 5 taajuuden siten, etta vastaanottimella 1 vastaanotetaan 
kulloinkin haluttuja taajuuksia. Vastaanotettavasta signaalista tallenne- 

35 taan naytteita naytevektoreiden muodostuslohkossa 12 eduliisesti 
kompleksisiksi naytevektoreiksi p(1), p(2)..,p(N). 
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Myos tassa suoritusmuodossa korrelointivaihe voidaan suorittaa osit- 
tain jo naytteenoton aikana, tai sen jalkeen kun N kappaletta naytevek- 
toreita p(l), p(2)...p(N) on muodostettu. Jos korrelointivaihe suoritetaan 
esim. siten, etta kunkin naytevektorin tallennuksen jalkeen lasketaan 
5 sille aika-taajuusmuunnos Hartley-muunnoksella, sopivimmin nopealla 
Hartley-muunnoksella (FHT), voidaan samaa aika-taajuusmuunninta 
kayttaa kaikissa N kappaleessa naytevektoreita p(1), p(2)..,p(N). Jos 
sen sijaan korrelointivaihe suoritetaan N naytevektorin tallennuksen jal- 
keen, on kaytettava joko kullekin naytevektorille omaa aika-taajuus- 

10 muunninta, tai aika-taajuusmuunnokset suoritetaan eri naytevektoreille 
perakkain samassa aika-taajuusmuuntimessa. Tassa suoritusmuodos- 
sa korrelointivaihe suoritetaan laskentapalvelimessa 26. Talloin elek- 
troniikkalaitteesta 24 lahetetaan edullisesti langattoman viestimen ra- 
dio-osan 22 avulla naytevektoreiden naytteiden arvot seka elektroniik- 

15 kalaitteen 24 aikatieto matkaviestinverkkoon. Matkaviestinverkossa 
nama tiedot vastaanotetaan tukiasemassa 25 ja valitetaan laskentapal- 
velimelle 26, jossa ne siirretaan kasiteltavaksi korrelointivaiheen suorit- 
tavaan lohkoon, jota kuvassa 6 esittaa tahdistuslohko 6. Naytevekto- 
reista p(1), p(2).„p(N) muodostetaan korrelaatiofunktiomatriisi C x , ku- 

20 ten keksinnon ensimmaisen edullisen suoritusmuodon mukaisen mene- 
telman kuvauksen yhteydessa on esitetty. 

Kullekin naytevektorille p(1), p(2)...p(N) suoritetaan diskreetti Hartley- 
muunnos 102, sopivimmin nopea Hartley-muunnos aikaisemmin tassa 
25 selityksessa esitetyn kaavan (4) mukaisesti. 

Laskentapalvelimessa 26 on tallennettuna sopivimmin kaikkien GPS- 
jarjestelman satelliittien (ei esitetty) C/A-koodia vastaavat vertailukoodit 
r(x), tai ne generoidaan laskentapalvelimessa 26. Korrelaatiovaiheessa 
30 valitaan tai generoidaan kulloinkin sen satelliitin vertailukoodi, jonka la- 
hettamaan signaaliin vastaanotinta tahdistetaan. Vertailukoodi on 
kaannetty ajallisesti takaperin. Talle takaperoiselle vertailukoodille, jota 
kuvassa 2 on merkitty viitteella r(x) , suoritetaan diskreetti Hartley- 



35 



muunnos 103, sopivimmin nopea Hartley-muunnos aikaisemmin tassa 
selityksessa esitetyn kaavan (5) mukaisesti. 
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Seuraavaksi korrelaatiovaiheessa suoritetaan kunkin naytevektorin p(i) 
Hartley-muunnostuloksen P(i) ja takaperoisen vertailukoodin F(jc) Hart- 

ley-muunnoksen R(x) valinen kertolasku aikaisemmin tassa selitykses- 
sa esitetyn kaavan (6) mukaisesti. 

5 

Naille kertolaskujen tuloksille suoritetaan viela kaanteinen Hartley- 
muunnos 105 (IFHT), jolloin tuloksena saadaan vertailukoodin r(x) ja 
vastaanotetun signaalin ristikorrelaatio kaikilla mahdollisilla kokonaislu- 
kuviiveilla (1023 kpl) (kaava 7). 

10 

Analyysivaiheessa korrelaatiofunktiomatriisin C x transpoosin kullekin 
riville suoritetaan Hartley-muunnos koherentin etsintamatriisin A x muo- 
dostamiseksi, jolloin voidaan suorittaa taajuusanalyysi todellisen taa- 
juussiirtyman selvittamiseksi. 

15 

Tahdistuslohkossa 6 voidaan tarvittaessa suorittaa viela epakoherentti 
summausvaihe signaalikohinasuhteen parantamiseksi. Epakoherentin 
summausvaiheen toteuttamiseksi toistetaan edella esitettyja naytevek- 
toreiden muodostusvaihetta, korrelointivaihetta ja analyysivaihetta K 
20 kertaa (ei esrtetty). Tama toistojen lukumaara K valitaan edullisesti si- 
ten, etta signaali-kohinasuhdetta saadaan riittavasti parannettua, kui- 
tenkin kohtuuliisessa ajassa. 

Selvitysvaiheessa (acquisition) tutkitaan etsintamatriisin ^alkioiden 
25 a x (/,y) arvoja ja pyritaan loytamaan 110 sellainen arvo, joka on toden- 

nakoisesti satelliitin lahettaman signaalin aiheuttama. Mikali tallainen 
arvo loytyy, laskentapalvelin 26 lahettaa tiedon koodivaiheesta ja taa- 
juuspoikkeamasta elektroniikkalaitteelle 24. Jos sen sijaan etsintamat- 
riisista A x el loydy tallaista arvoa, eli tutkitulla taajuusalueella ei toden- 

30 nakoisesti vastaanotettu haettavan satelliitin lahettamaa signaalia, 
muutetaan tutkittavaa taajuuskaistaa ja suoritetaan edella esitetyt vai- 
heet etsintamatriisin muodostamlseksi. Taman aikaansaamiseksi las- 
kentapalvelin 26 lahettaa elektroniikkalaitteelle 24 tiedon siita, etta tut- 
kittavana olleella taajuusalueella ei havaittu satelliitin signaalia, jolloin 

35 elektroniikkalaitteen vastaanotin 1 sopivimmin muuttaa tutkittavaa taa- 
juusaluetta ja aloittaa edella esitettyjen vaiheiden suorittamisen talla 
muutetulla taajuusalueella. Talla menetelmalla saadaan koko tutkittava 
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6 kHz:n alue kaytya lapi toistamalla edella esitetyt vaiheet tarvittavan 
rnonta kertaa. 

Tarvittaessa voidaan edella esitetyt vaiheet toistaa koko tutkittavalle 
5 taajuuskaistalle ja tallentaa laskentapah/elimessa 26 eri toistokerroilla 
rnuodostetut etsintamatriisit tai vain mahdolliset huippukohdat, ennen 
kuin etsitaan suurinta korrelaatiohuippua. Talloin voidaan pienentaa 
virhetulkintojen mahdollisuutta mm. sellalsessa tilanteessa, jossa kyn- 
nysarvo on asetettu liian pieneksi ja hairiosignaali voi arheuttaa virhe- 
10 tulkinnan. 

Tama keksinnon eraan toisen edullisen suoritusmuodon mukainen me- 
netelrna on erityisen tehokas mm. siita syysta, etta sijainninmaarityk- 
sessa kaytettavan laskentapalvelimen 26 sijaintimaarityksissa on erit- 
15 tain pieni epatarkkuus, joka kaytannossa on alle millimetrin suuruus- 
luokkaa. Tama johtuu siita, etta sijainninmaSrityspalvelin sijaitsee tar- 
kasti tunnetussa paikassa (esim. tukiasema) ja silla on tavallisesti kay- 
tettavissa satelliittien sen hetkinen efemeris-tieto ja almanakkatieto, 

20 Verkkopohjaisen palkannuksen eras perusidea on sils toteuttaa kaikki 
vaativa laskenta tiedonsiirtoverkkoon sijoitetussa laskentapalvelimes- 
sa 26 tai vastaavassa laitteessa. Elektroniikkalaitteen vastaanottimen 1 
tehtavana on talloin lukittua vahintaan neljan satelllitin signaaliin ja lilt— 
taa niiden koodisignaaleihin aikatieto pseudoetaisyyksien laskemista 

25 varten. Taman jalkeen elektroniikkalaite 24 IShettaa signaalit laskenta- 
palvelimelle 26 varsinaisen sijainnin laskemiseksi. 

Verkkopohjaisessa paikannuksessa vastaanotin 1 koostuu minimissaan 
vain tarvittavasta laitteistosta alibitti- ja koodivaihemittausten tekemi- 
30 seen ja niiden valittamiseen verkon kautta palvelimelle. 

Vaikka edella keksintda on selostettu paaasiassa vastaanottimen tah- 
distumisen yhteydessa, on selvaa, etta keksintoa voidaan soveltaa 
myos tahdistumisen yllapitamisen yhteydessa. 

35 

Nyt esilla olevaa keksintoa ei ole rajoitettu ainoastaan edella esitettyihin 
suoritusmuotoihin, vaan sita voidaan muunnella oheisten patenttivaati- 
musten puitteissa. 
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Patenttivaatimukset: 

1. Menetelma vastaanottimen (1) tahdistustoimintojen vastaanotettuun 
koodimoduloituun hajaspektrisignaaliin suorittamiseksi, jossa menetel- 

5 mSssa kaytetaan ainakin yhta vertailukoodia (r(x)), joka vastaa jotakin 
moduloinnissa kaytettya koodia, selvitetaan vastaanotetun signaalin 
taajuussiirtyma ja moduloinnissa kaytetyn koodin koodivaihe, tunnettu 
siita, etta menetelmassa suoritetaan ainakin seuraavia vaiheita: 

- naytevektoreiden muodostusvaihe, jossa vastaanotetusta signaalista 
10 otetaan naytteita (101) naytevektoreiden (p(1), p(2)...p(N)) muodos- 

tamiseksi, 

- korrelaatiovaihe(102, 103, 104, 105), jossa muodostetaan mainitun 
vertailukoodin (r(x)) perusteella ensimmainen Hartley-muun- 
nos (R(x)), ja kunkin naytevektorin (p(1), p(2),..p(N)) perusteella toi- 

15 nen Hartley-muunnos (P(i)), suoritetaan mainitun vertailukoodin pe- 
rusteella muodostetun ensimmaisen Hartley-muunnoksen (R(x)) ja 

kunkin naytevektorin (p(1), p(2)„.p(N)) perusteella muodostetun toi- 
sen Hartley-muunnoksen (P(i)) valinen kertolasku, suoritetaan kulle- 
kin kertolaskutulokselle (M x (i)) kaanteinen Hartley-muunnos, ja 
20 - selvitysvaiheCHO), jossa selvitetaan kertolaskutulosten (M x (i)) 
kaanteisten Hartley-muunnosten (m x (i)) perusteella taajuussiirtyma ja 
koodivaihe. 

2. Patenttivaatimuksen 1 mukainen menetelma, tunnettu siita, etta 
25 menetelmassa maaritetaan tutkittava taajuusalue, jaetaan maaritetty 

taajuusalue kahteen tai useampaan osaan, jolloin kullekrn osalle suori- 
tetaan mainitut naytevektoreiden muodostusvaihe ja korrelaatiovaihe, ja 
etta menetelmassa suoritetaan lisaksi analyysivaihe, jossa tallennetaan 
kertolaskutulosten (M x (i)) kaanteisten Hartley-muunnosten (m x (i)) arvoja 
30 koherentin etsintamatriisin (A x ) muodostamiseksi, jolloin selvrtysvaihe 

suoritetaan maaritetyn taajuusalueen tutkimisen jalkeen, ja taajuussiir- 
tyma ja koodivaihe maaritetaan koherentin etsintamatriisin (A x ) suu- 

rimman arvon perusteella. 
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3. Patenttivaatimuksen 2 mukainen menetelma, tunnettu siita, etta 
menetelmassa maaritetaan kynnysarvo, jolloin selvitysvaiheessa kayte- 
taan taajuussiirtyman ja koodivaiheen selvittamiseksi mainltun kynnys- 

5 arvon ylittavia koherentin etsintamatriisin (A x ) alkioiden (a x (/,/)) suu- 

ruusarvoja. 

4. Patenttivaatimuksen 2 tai 3 mukainen menetelma, tunnettu siita, 
etta menetelmassa toistetaan (1 08) mainittuja naytevektoreiden muo- 

10 dostusvaihetta, korrelaatiovaihetta ja analyysivaihetta ainakin kahden 
koherentin etsintamatriisin (A x ) muodostamiseksi, ja etta menetel- 
massa lisaksi suoritetaan summausvaihe (109), jossa muodostetaan 
epakoherentti etsintamatriisi (s x ) summaamalla epakoherentisti kulla- 
kin toistokerralla muodostetut koherentin etsintamatriisin (A x ) vastinal- 

15 kioiden arvot (a x (/,/)), jolloin mainitussa selvitysvaiheessa kaytetaan 
mainittua epakoherenttia etsintamatriisia (s x ) taajuussiirtyman ja koodi- 
vaiheen maarityksessa. 

5. Patenttivaatimuksen 4 mukainen menetelma, tunnettu siita, etta 
20 menetelmassa maaritetaan tutkittava taajuusalue, jaetaan maaritetty 

taajuusalue kahteen tai useampaan osaan, jolloin kullekin osalle suori- 
tetaan mainitut naytevektoreiden muodostusvaihe, korrelaatiovaihe, 
analyysivaihe, ja summausvaihe, seka tallennetaan epakoherentin et- 
sintamatriisin (s x ) alkioiden arvoja, jolloin selvitysvaihe suoritetaan 

25 maaritetyn taajuusalueen tutkimisen jalkeen, ja taajuussiirtyma ja koo- 
divaihe maaritetaan suurimman arvon perusteella, 

6- Patenttivaatimuksen 4 tai 5 mukainen menetelma, tunnettu siita, 
etta menetelmassa maaritetaan kynnysarvo, jolloin selvitysvaiheessa 
30 kaytetaan taajuussiirtyman ja koodivaiheen selvittamiseksi mainitun 
kynnysarvon ylittavia epakoherentin etsintamatriisin (s x ) alkioiden suu- 
ruusarvoja, 

7. Jonkln patenttivaatimuksen 1 — 6 mukainen menetelma, tunnettu 
35 siita, etta korrelaatiovaiheessa kaytetaan mainittua vertailukoodia (r(x)) 
vastaavaa takaperoista koodia (r(x)) ensimmaisen Hartley-muunnok- 
sen (R(x)) muodostamiseksi. 
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8. Jonkin patenttivaatirnuksen 1 — 6 mukainen menetelma, tunnettu 
siita, etta korrelaatiovaiheessa kdytetaan kutakin naytevekto- 
ria (p(1), p(2)„.p(N)) vastaavaa takaperotsta koodla toisen Hartley- 

5 muunnoksen (P(/)) muodostamiseksi. 

9. Paikannusjarjestelma, joka kasittaa ainakin vastaanottimen (1), tah- 
distusvalineet (6) vastaanottimen (1) tahdistustormintojen suorittarni- 
seksi vastaanotettuun koodimoduloituun hajaspektrisignaaliin, valineet 

10 (16) ainakin yhden vertailukcrodin (r(x)) kayttarniseksi tahdistamisen 
yhteydessa, joka vertailukoodi (r(x)) vastaa jotakin moduloinnissa kay- 
tettya koodia, ja valineet (15) vastaanotetun signaalin taajuussiirtyman 
ja moduloinnissa kaytetyn koodin koodivaiheen selvittamiseksi, tun- 
nettu siita, etta tahdistusvalineet kasittavat ainakin: 

15 - naytevektoreiden muodostusvalineet (12, 16) naytevektorei- 
den (p(1), p(2)...p(N)) muodostamiseksi vastaanotetusta signaalista, 

- korrelaatiovalineet (6, FHT1, FHT2...FHTN) mainitun vertailukoo- 
din (r(x)) perusteella ensimmaisen Hartley-muunnoksen (R(x)) muo- 
dostamiseksi, ja kunkin naytevektorin (p(1), p(2).,.p(N)) perusteella 

20 toisen Hartley-muunnoksen (P(i)) muodostamiseksi, valineet maini- 
tun vertailukoodin perusteella muodostetun ensimmaisen Hartley- 
muunnoksen (R(x)) ja kunkin naytevektorin (p(1) f p(2)...p(N)) perus- 
teella muodostetun toisen Hartley-muunnoksen (P(i)) valisen kerto- 
laskun suorittamiseksi, ja valineet kullekin kertolaskutulokselle (M x (i)) 

25 kaanteisen Hartley-muunnoksen suorittamiseksi korrelaatiofunktio- 
matriisin (C x ) muodostamiseksi, 

- selvitysvalineet (7, 17) taajuussiirtyman ja koodivaiheen selvittami- 
seksi kayttarnalla kertolaskutulosten (M x (i)) kaanteisten Hartley- 
muunnosten (m,c(i)) arvoja. 

30 

10. Patenttivaatirnuksen 9 mukainen paikannusjarjestelma, tunnettu 
siita, etta vastaanotetut signaalit ovat GPS-jarjestelman satelliittien la- 
hettamia signaaleja. 
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11. Patenttivaatimuksen 9 tai 10 mukainen paikannusjarjestelma, joka 
kasittaa tiedonsiirtoverkon, tunnettu siita, etta ainakin osa tahdistus- 
valineista (6) on muodostettu tiedonsiirtoverkon yhteyteen, jolloin tie- 
donsiirtoverkon ja vastaanottimen (1) valille on jarjestetty muodostetta- 
vaksi tiedonsiirtoyhteys. 

12. Patenttivaatimuksen 11 mukainen paikannusjarjestelma, tunnettu 
siita, etta tiedonsiirtoverkko kasittaa ainakin matkaviestinverkon. 

1 3. Patenttivaatimuksen 9 tai 1 0 mukainen paikannusjarjestelma, tun- 
nettu siita, etta tahdistusvalineet (6) on muodostettu vastaanottimeen 
0). 

15 14. Vastaanotin (1), joka kasittaa ainakin tahdistusvalineet (6) vas- 
taanottimen (1) tahdistustoimintojen vastaanotettuun koodimoduloituun 
hajaspektrisignaaliin suorittamiseksi, jossa vastaanottimessa (1) on va- 
lineet (16) ainakin yhden vertailukoodin (r(x)) kayttamiseksi tahdistami- 
sen yhteydessa, joka vertailukoodi (r(x)) vastaa jotakin moduloinnissa 
20 kaytettya koodia, ja valineet (15) vastaanotetun signaalin taajuussiirty- 
man ja moduloinnissa kaytetyn koodin koodivaiheen selvlttamiseksi, 
tunnettu siita, etta tahdistusvalineet kasittavat ainakin: 
- naytevektoreiden muodostusvalineet (12, 16) naytevektorei- 
den (p(l), p(2)...p(N)) muodostamiseksi vastaanotetusta signaalista, 
25 - korrelaatiovalineet (6, FHT1 , FHT2...FHTN) malnitun vertailukoo- 
din (r(x)) perusteella ensimmaisen Hartley-muunnoksen (R(x)) muo- 
dostamiseksi, ja kunkin naytevektorin (p(1), p(2)...p(N)) perusteella 
toisen Hartley-muunnoksen (P(i)) muodostamiseksi, valineet maini- 
tun vertailukoodin perusteella muodostetun ensimmaisen Hartley- 
30 muunnoksen (R(*)) ja kunkin naytevektorin (p(1), p(2)...p(N)) perus- 
teella muodostetun toisen Hartley-muunnoksen (P(i)) valisen kerto- 
laskun suorittamiseksi, ja valineet kullekin kertolaskutulokselle (M x (i)) 
kaanteisen Hartley-muunnoksen suorittamiseksi korrelaatiofunktio- 
matriisin (C x ) muodostamiseksi, 

35 - selvitysvalineet (7, 17) taajuussilrtyman ja koodivaiheen selvittami- 
seksi kayttamalla kertolaskutulosten (M x (i)) kaanteisten Hartley- 
muunnosten (mx(i)) arvoja. 
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15. Patenttivaatimuksen 14 mukainen vastaanotin (1), funneiffiu siita, 
etta se kasittaa lisaksi ainakin 

- valineet (17) tutkittavan taajuusalueen maarittamiseksi, valineet (5 f 8, 

9) maaritetyn taajuusalueen jakamiseksi kahteen tai useampaan 

5 osaan, jolloin kullekin osalle on jarjestetty suoritettavaksi mainitut 
naytevektoreiden muodostus ja korrelaatiofunktiomatriisin (C x ) muo- 

dostus, 

- valineet koherentin etsintamatriisin (A x ) muodostamiseksi, 

- valineet (16) koherentin etsintamatriisin (A x ) alkioiden (a Jf (/\y)) arvo- 
10 jen tallentamiseksi, ja 

- valineet (15) taajuussiirtyman ja koodivaiheen maarittamiseksi kohe- 
rentin etsintamatriisin (A x ) suurimman arvon perusteella. 

16. Patenttivaatimuksen 14 tai 15 mukainen vastaanotin (1), fyironetlty . 
15 siita, etta se kasittaa valineet (15) kynnysarvon maarittamiseksi, ja vali- 
neet (15) mainitun kynnysarvon ja mainittujen koherentin etsintamatrii- 
sin (A x ) arvojen vertaamiseksi taajuussiirtyman ja koodivaiheen selvit- 

tamiseksi. 

20 17. Patenttivaatimuksen 14, 15 tai 16 mukainen vastaanotin (1), tun- 
nettu siita, etta mainittuja naytevektoreiden muodostusta, korrelaatio- 
funktiomatriisin (C x ) muodostusta ja koherentin etsintamatriisin (A x ) 

muodostusta on jarjestetty toistettavaksi koherentin etsintamatriisin 
(A Kk ) muodostamiseksi ainakin kaksi kertaa p ja etta vastaanotin (1) 

25 kasittaa lisaksi summausvalineet epakoherentin etsintamatriisin (s x ) 

muodostamiseksi summaamalla kullakin toistokerralla muodostetut ko- 
herentin etsintamatriisin (A x ) vastinalkioiden arvot (a Xr ^(/ r y)) I jolloin 

mainitussa selvrtysvaiheessa kaytetaan mainittua epakoherenttia etsin- 
tamatriisia (S x ) taajuussiirtyman ja koodivaiheen maarityteessa. 

30 

18. Patenttivaatimuksen 14 mukainen vastaanotin (1). funnetty siita, 
etta se kasittaa lisaksi ainakin 

- valineet (17) tutkittavan taajuusalueen maarittamiseksi, 

- valineet (5, 8, 9) maaritetyn taajuusalueen jakamiseksi kahteen tai 

35 useampaan osaan, jolloin kullekin osalle on jarjestetty suoritettavaksi 

mainitut naytevektoreiden muodostus ja korrelaatiofunktiomatrii- 
sin (C x ) muodostus, 
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- valineet koherentin etsintamatriisin (A x ) muodostamiseksi, 

- valineet koherentin etsintamatriisin (A x ) summaamiseksi epakohe- 
renttiin etsintamatriisiin (5 X ), 

- valineet (16) epakoherentin etsintamatriisin (s x ) alkioiden arvojen 
5 tallentamiseksi, ja 

- valineet (15) taajuussiirtyman ja koodivaiheen maarittamiseksi suu- 
rimman arvon perusteella. 

19. Patenttivaatimuksen 18 mukainen vastaanotin (1), Himirora^y siita, 
10 etta se kasittaa valineet (15) kynnysarvon maarittamiseksi, ja valineet 
(15) mainitun kynnysarvon ja mainitun epakoherentin etsintamatriisin 
(S x ) arvojen vertaamiseksi taajuussiirtyman ja koodivaiheen selvittami- 
seksi. 

15 20. Jonkin patenttivaatimuksen 14 — 19 mukainen vastaanotin (1), 
twranetty siita, etta korrelaatiovalineet kasittavat valineet (16) mainfttua 
vertailukoodia (r(x)) vastaavan takaperoisen koodin (f(x)) Hartley- 
muunnoksen (R(x)) muodostamiseksi- 

20 21. Jonkin patenttivaatimuksen 14 — 19 mukainen vastaanotin (1), 
turonettu siita, etta korrelaatiovalineet kasittavat valineet (16) kutakin 
naytevektoria (p(l), p(2).,.p(N)) vastaavan takaperoisen koodin Hartley- 
muunnoksen (P(/)) muodostamiseksi. 

25 22. Elektroniikkalaite (24), joka kasittaa ainakin sijainninmaaritysvas- 
taanottimen (1), tahdistusvalineet (6) sijainninmaaritysvastaanotti- 
men (1) tahdistustoimintojen lahetettyyn koodimoduloituun hajaspektri- 
signaaliin suorittamiseksi, jossa sijainninmaaritysvastaanottimessa (1) 
on valineet (16) ainakin yhden vertailukoodin (r(x)) kayttamiseksi 

30 tahdistamisen yhteydessa, joka vertailukoodi (r(x)) vastaa jotakin modu- 
loinnissa kaytettya koodia, ja joka elektroniikkalarte kasittaa valineet 
(15) lahetetyn signaalin taajuussiirtyman ja moduloinnissa kaytetyn 
koodin koodivaiheen selvittamiseksi, fuiminiettiui siita, etta tahdistusvali- 
neet kasittavat ainakin: 

35 - naytevektoreiden muodostusvSlineet (12, 16) naytevektorei- 
den (p(1), p(2)„.p(N)) muodostamiseksi vastaanotetusta signaalista, 
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- korrelaatiovalineet (6, FHT1 , FHT2...FHTN) mainitun vertailukoo- 
din (r(x)) perusteella ensimmaisen Hartley-muunnoksen (R(x)) muo- 

dostamiseksi, ja kunkin naytevektorin (p(l), p(2)...p(N)) perusteella 
toisen Hartley-muunnoksen (P(i)) muodostamiseksi, valineet maini- 
5 tun vertailukoodin perusteella muodostetun ensimmaisen Hartley- 
muunnoksen (R(x)) ja kunkin naytevektorin (p(l), p(2)...p(N)) perus- 
teella muodostetun toisen Hartley-muunnoksen (P(i)) valisen kerto- 
laskun suorittamiseksi, ja valineet kullekin kertolaskutulokselle (M x (i)) 

kaanteisen Hartley-muunnoksen suorittamiseksi korrelaatiofunktio- 
10 matriisin (C x ) muodostamiseksi, ja 

- selvitysvalineet (7, 17) taajuussiirtyman ja koodivaiheen selvittami- 
seksi kayttamalla kertolaskutulosten (M s (i)) kaanteisten Hartley- 
muunnosten (m x (i))arvoja. 

15 23. Patenttivaatimuksen 22 mukainen elektroniikkalaite (24), joka kasit- 
taa valineet (14, 15, 17) elektroniikkalaitteen (24) sijainnin selvittami- 
seksi ja valineet (18) sijaintitiedon tallentamiseksi, funnaffitu siita, etta 
elektroniikkalaite (24) kasittaa lisaksi valineet (17) tutkittavan taajuus- 
alueen maarittamiseksi, ja valineet (1 7, 22, 23) aloitustaajuuden valit- 

20 semiseksi mainitulta taajuusalueelta vastaanottimeen (1) tallennetun 
sijaintitiedon perusteella. 

24. Patenttivaatimuksen 22 tai 23 mukainen elektroniikkalaite (24), 
funnettu siita, etta korrelaatiovalineet kasittavat valineet (16) mainittua 

25 vertallukoodia (r(x)) vastaavan takaperoisen koodin (f(x)) Hartley- 
muunnoksen (R(*)) muodostamiseksi. 

25. Patenttivaatimuksen 22 tai 23 mukainen elektroniikkalaite (24), 
tyiramieffiu siita, etta korrelaatiovalineet kasittavat valineet (16) kutakin 

30 naytevektoria (p(1 ), p(2)...p(N)) vastaavan takaperoisen koodin Hartley- 
muunnoksen (P(/)) muodostamiseksi. 

26. Jonkin patenttivaatimuksen 22—25 mukainen elektroniikka- 
laite (24), tunnettu siita, etta se kasittaa lisaksi valineet tiedonsiirto- 

35 toimintojen suorittamiseksi. 
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27. Patenttivaatimuksen 26 mukainen elektroniikkalaite (24), tunnettu 
siita, etta se kasittaa valineet (22, 23) tiedonsiirtoyhteyden muodosta- 
miseksi tiedonsiirtoverkkoon, jolloin valineet elektroniikkalaltteen (24) 
sijainnin servittamiseksl kasittavat valineet (22, 23) sijainninmaarityk- 
5 sessa tarvittavan tiedon lahettamiseksi tiedonsiirtoverkkoon, ja vali- 
neet (22, 23) sijainninmaarityksessa kaytettavan tiedon hakemiseksi 
tiedonsiirtoverkosta, jolloin ainakin osa sijainninrnaaritystoiminnoista on 
jarjestetty suoritettavaksi tiedonsiirtoverkossa. 

10 28. Patenttivaatimuksen 26 tai 27 mukainen elektroniikkalaite (24), 
tunnettu siita, etta valineet tiedonsiirtotoimintojen suorittamiseksi kasit- 
tavat ainakin valineet matkaviestintoimintojen suorittamiseksi. 

29. Elektroniikkalaite (24), joka kasittaa ainakin vastaanottimen (1), ja 
15 joka elektroniikkalaite (24) on tarkoitettu kaytettavaksi paikannusjarjes- 
telman yhteydessa, joka kasittaa ainakin tahdistusvalineet (6) vas- 
taanottimen (1) lahetettyyn koodimoduloituun hajaspektrisignaaiiin tah- 
distustoimintojen suorittamiseksi, valineet (16) ainakin yhden vertailu- 
koodin (r(x)) kayttamiseksi tahdistamisen yhteydessa, joka vertailukoodi 
20 (r(x)) vastaa jotakin moduloinnissa kaytettya koodia, valineet (15) lahe- 
tetyn signaalin taajuussiirtyman ja moduloinnissa kaytetyn koodin koo- 
divaiheen selvittamiseksi, ja tiedonsiirtoverkon, tunnettu siita, etta 
elektroniikkalaite (24) kasittaa lisaksi ainakin naytevektoreiden muodos- 
tusvalineet (12, 16) naytevektoreiden (p(1), p(2)...p(N)) muodostami- 
25 seksi vastaanotetusta signaalista, ja lahetysvalineet (22, 23) muodos- 
tettujen naytevektoreiden (p(1), p(2)...p(N)) ja aikatiedon lahettamiseksi 
tiedonsiirtoverkkoon, jolloin paikannusjarjestelma kasittaa lisaksi aina- 
kin 

- korrelaatiovalineet (6, FHT1, FHT2...FHTN) mainitun vertailukoo- 
30 din (r(x)) perusteella ensimmaisen Hartley-muunnoksen (R(x)) muo- 

dostamiseksi, ja kunkin naytevektorin (p(1), p(2)...p(N)) perusteella 
toisen Hartley-muunnoksen (P(i)) muodostamiseksi, valineet maini- 
tun vertailukoodin perusteella muodostetun ensimmaisen Hartley- 
muunnoksen (R(x)) ja kunkin naytevektorin (p(1), p(2)...p(N)) perus- 

35 teella muodostetun toisen Hartley-muunnoksen (P(i)) valisen kerto- 
laskun suorittamiseksi, ja valineet kullekin kertolaskutulokselle (M x (i)) 
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kaanteisen Hartley-muunnoksen suorittamiseksi korrelaatiofunktio- 
matriisin (C x ) muodostamiseksi, ja 

- selvitysvaiineet (7, 17) taajuusslirtyman ja koodivaiheen selvittami- 
seksi kayttamalla kertolaskutulosten (M x (i)) kaanteisten Hartley- 
5 muunnosten (m x (i))arvoja. 

30. Patenttivaatimuksen 29 imukainen elektroniikkalaite (24), Hyiiniinietttty 
siita, etta se kasittaa lisaksi vastaanottovalineet (23, 22) tiedon selvlte- 
tysta taajuussiirtymasta ja koodivaiheesta vastaanottamiseksi tiedonsiir- 

10 toverkosta. 

31 . Patenttivaatimuksen 29 tai 30 mukainen elektroniikkalaite (24), joka 
on tarkoitettu kaytettavaksi paikannusjarjestelman yhteydessa, joka 
kasittaa lisaksi ainakin valineet (26) elektroniikkalaitteen sijainnin maa- 

15 rittamiseksi, ftymnmefiSy siita, etta elektroniikkalaite (24) kasittaa lisaksi 
valineet (22, 23) sijaintitiedon vastaanottamiseksi tiedonsiirtoverkosta. 

32. Patenttivaatimuksen 29, 30 tai 31 mukainen elektroniikkalaite (24), 
tuoimineitty siita, etta valineet tiedonsiirtotoimintojen suorittamiseksi kasit- 

20 tavat ainakin valineet matkaviestintoimintojen suorittamiseksi. 
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(57) Tiivistelma : 

Keksinnon kohteena on menetelma vastaanottimen (1) 
tahdistustoimintojen vastaanotettuun koodimoduloituun 
hajaspektrisignaaliin suorittamiseksi. ' Menetelmassa 
kaytetaan ainakin yhta vertailukoodia (r(x)), joka vaslaa 
jot akin moduloinnissa kaytettya koodia. Menetelmassa 
selvitetaan vastaanotetun signaalin taajuussiirtyma ja 
moduloinnissa kaytetyn koodin koodivaihe. Menetelmas- 
sa suoritetaan lisaksi ainakin seuraavia vaiheita: 

- naytevektoreiden muodostusvaihe, jossa vastaanote- 
tusta signaalista otetaan naytte'rta (101) naytevektorei- 
den (p(1), p(2)...p(N)) muodostamiseksi, 

- korrelaatiovaihe (102, 103, 104, 105), jossa muodos- 
tetaan mainitun vertailukoodin (r(x)) perusteella en- 
simmainen Hartley-muunnos (R(x)), ja kunkin nayte- 

vektorin (p(1), p(2)...p(N)) perusteella toinen Hartley- 
muunnos (P(i)), suoritetaan mainitun vertailukoodin 
perusteella muodostetun ensimmaisen Hartley-muun- 
noksen (R(x)) ja kunkin naytevekto- 

rin (p(1), p(2)...p(N)) perusteella muodostetun toisen 
Hartley-muunnoksen (P(i)) valinen kertolasku, suorite- 
taan kullekin kertolaskutulokselle (M x (i)) kaanteinen 
Hartley-muunnos, ja 

- selvitysvaihe(HO), jossa selvitetaan kertolaskutulos- 
ten (Mx(i)) kaanteisten Hartley-muunnosten (m,t(i)) pe- 
rusteella taajuussiirtyma ja koodivaihe. 

Fig.1 
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Method for synchronizing a receiver, a positioning system, a receiver and an 
electronic device 

The present invention relates to a method according to the preamble of 
Claim 1 for performing the synchronization operations of a receiver, a 
positioning system according to the preamble of Claim 9, a receiver 
according to the preamble of Claim 14, an electronic device according to the 
preamble of Claim 22, and an electronic device according to the preamble of 
Claim 29. 

A well-known positioning system is the Global Positioning System, GPS, 
which now comprises over 20 satellites, 4 of which or more are 
simultaneously visible to a receiver. In Finland, for instance, it is possible to 
detect even more than 14 satellites simultaneously, depending on the 
latitude, due to the visibility over the North Pole. These satellites transmit, 
for instance, Ephemeris data and data of the satellite time. Normally, the 
receiver used in positioning determines its location by calculating the transit 
time to the receiver of a signal transmitted simultaneously from several 
satellites of the positioning system. For determining its location, the receiver 
must typically receive the signal of at least four visible satellites to be able to 
calculate the location. 

Each satellite of the GPS system in operation transmits a so-called L1 signal 
at the carrier frequency of 1575.42 MHz. This frequency is also denoted 
154f 0 , where f 0 =10.23 MHz. In addition, the satellites transmit an L2 signal 
at the carrier frequency of 1227.6 MHz, or 120f 0 . In the satellite, these 
signals are modulated with at least one pseudo-random sequence. Each 
satellite has a different pseudo-random sequence. After the modulation, the 
signal is a code modulated wideband signal. The modulation technique 
used makes it possible for the receiver to distinguish between signals 
transmitted by different satellites, although the carrier frequencies used in 
the transmission are essentially the same. This modulation technique is 
called Code Division Multiple Access (CDMA). The pseudo-random 
sequence used in each satellite for the modulation of the L1 signal is, for 
instance, a so-called C/A code (Coarse/Acquisition code), an example of 
which is the Gold code. Each GPS satellite transmits a signal by using a 
unique C/A code. The codes are formed as the modulo 2 sum of two 1023- 
bit binary sequences. The first binary sequence G1 has been formed with 
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the polynome X 10 +X 3 +1, and the second binary sequence G2 has been 
formed by delaying the polynome X 10 +X 9 +X 8 +X 6 +X 3 +X 2 +1 so that each 
satellite has a different delay. This arrangement enables generating different 
C/A codes with a similar code generator. The C/A codes are binary codes, 
5 the chipping rate of which in the GPS system is 1 .023 MHz. The C/A code 
comprises 1023 chips, whereby the repeating time of the code, epoch, is 1 
ms. The carrier wave of the L1 signal is further modulated with navigation 
information at a bit rate of 50 bit/s. The navigation information comprises 
information about the "health", orbit and time of the satellite, for instance. 

10 

The satellites monitor the condition of their equipment during operation. For 
example, the satellites can use so-called watch-dog functions to detect 
errors in the equipment and notify about them Errors and malfunctions can 
be either momentary or long-term in nature. On the basis of the health 

15 information, some of the errors can be possibly compensated for, or the 
information transmitted by a malfunctioning satellite can be entirely ignored. 
In addition, in a situation where the signal of more than four satellites can be 
received, the information received from different satellites can be 
emphasized in different ways on the basis of the health information. Thus it 

20 is possible to minimize errors that satellites which seem to be unreliable 
may cause in the measurements. 

In order to detect the signals of the satellites and to identify the satellites, 
the receiver must perform synchronization, whereby the receiver seeks out 
25 the signal of each satellite each time and tries to synchronize with this signal 
so that the data transmitted with the signal can be received and 
demodulated. 

A location determination receiver must carry out synchronization when, for 
30 instance, the receiver is switched on, and also in a situation where the 
receiver has not been able to receive the signal of any satellite for a long 
time. In portable devices, for example, a situation like this can arise easily, 
because the device is moving, and the antenna of the device is not always 
in an optimal position in relation to the satellites, which weakens the 
35 strength of the signal coming to the receiver. In city areas, buildings have a 
blocking effect on the received signal, and in addition, the signal can 
proceed to the receiver along several different routes, such as straight from 
the satellite (line-of-sight) and as reflected from buildings. As a result of this 
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multipath propagation, the same signal is received as several signals in 
different phases. 

The positioning arrangement has two main functions: 

5 

1. calculating the pseudorange of the receiver to different GPS 
satellites, and 

2. determining the location of the receiver using the calculated 
pseudoranges and the location information of the satellites. The 

10 location information of the satellites at any given moment can be 

calculated on the basis of the Ephemeris and time correction 
information received from the satellites. 

The distances to the satellites are called pseudoranges, because the time is 
15 not known accurately in the receiver. Because the time is not known with an 
absolute accuracy, the location and time must be determined preferably by 
iterating the measured data with a linearized equation group. Then the 
determination of the location and time is repeated, until a sufficient accuracy 
of time and location has been achieved. 

20 

The calculation of the pseudorange can be carried out by measuring the 
mutual, apparent transit time delays of the signals of different satellites. 
After the receiver has become synchronized with the received signal, the 
information transmitted in the signal is detected. 

25 

Almost all known GPS location determination receivers (GPS navigation 
receivers) use correlation methods for carrying out synchronization 
functions,. i.e. code acquisition and/or maintaining synchronization (tracking). 
The reference codes r(x) or pseudo random sequences of different satellites 

30 have been saved or are generated locally in the location determination 
receiver. A down conversion is performed on the received signal, whereafter 
the receiver multiplies the received signal with the saved pseudo random 
sequence. The signal formed as the result of the multiplication is integrated 
or low-pass filtered, whereby the result tells if the received signal contained 

35 a signal transmitted by the satellite. The multiplication carried out in the 
receiver is repeated so that each time the phase of the pseudo random 
sequence saved in the receiver is shifted. The right phase is concluded from 
the correlation result preferably so that when the correlation result is the 
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highest, the right phase has been found. Then the receiver is correctly 
synchronized with the received signal. 

After the synchronization with the code has been performed, the frequency 
5 is adjusted and phase locking is carried out. This correlation result also 
expresses the information transmitted in the GPS signal. 

The above mentioned synchronization and frequency adjustment process 
must be repeated on each signal of a satellite which is received in the 

10 receiver. Therefore, this process consumes a lot of time especially in a 
situation where the received signals are weak. In order to make this process 
faster, several correlators are used in some prior art receivers, whereby 
more correlation peaks can be searched for simultaneously. Only by 
increasing the number of correlators it is not possible to make the 

15 synchronization and frequency adjustment process much faster in practical 
applications, because the number of correlators cannot be increased to 
infinity. 

Some prior art GPS receivers use the FFT technique for determining the 
20 Doppler shift of the received GPS signal in connection with ordinary 
correlators. In these receivers, correlation is used to reduce the bandwidth 
of the received signal to 10kHz-30kHz. This narrow-band signal is 
analyzed with FFT algorithms to determine the carrier frequency. 

25 A GPS receiver and a method for processing GPS signals have been 
presented in the international patent application WO 97/14057. In the main, 
the receiver presented in this specification comprises two separate 
receivers, of which the first receiver is intended for use in a situation where 
the strength of the received signal is sufficient, and the second receiver is 

30 intended for use in a situation where the strength of the received signal is 
not sufficient for a sufficiently accurate location determination when the first 
receiver is used. In this second receiver, the received signal is digitized and 
saved in memory media, whereby these saved signals are later processed 
in a digital signal processing unit. The digital signal processing unit performs 

35 convolution operations to the received, digitized signal. The purpose of 
these convolution operations is to calculate pseudoranges. The number of 
PM frames saved in the memory media is typically from one hundred to one 
thousand, which corresponds to a signal of the length of 1 00 ms to 1 s. After 



this, a saved code corresponding to the code of the satellite to be examined 
is retrieved from the memory of the receiver to be used for analyzing the 
received signal. 

The Doppler shift is also removed in the receiver. The amount of this 
Doppler shift has been determined either with the first receiver or on the 
basis of information received from a base station of the GPS system. After 
this, a coherent summing of consecutive frames is performed. A fast Fourier 
transformation is performed on the data set received as the result of this 
summing. A multiplication is performed on the Fourier transformation result, 
using the complex conjugate of the Fourier transformation of the reference 
signal saved in the memory media as one multiplier. An inverse Fourier 
transformation is further performed on the result of this multiplication, 
whereby a number of correlation results are received. In this specification, 
correlation has been replaced by a Fourier transformation, whereby the 
number of calculations has been reduced. According to the specification, the 
method makes location determination 10 to 100 times faster as compared to 
the solutions known at the time of filing the specification. 

Matched filters are used in some prior art receivers for the locking and 
tracking of the signal. However, the use of matched filters entails the 
drawback that they consume a relatively high amount of power, which limits 
the use of matched filters especially in mobile devices. 

It is an objective of the present invention to provide a receiver, in which 
synchronization with the transmitted signal can be performed substantially 
faster than in the prior art receivers. The invention is especially suitable for 
use in location determination receivers, but also in other receivers, 
preferably in CDMA receivers, in which the receiver must become 
synchronized with a spread spectrum signal. The invention is based on the 
idea that the time-to-frequency conversions and the inverse or frequency-to- 
time conversions are performed with a fast Hartley conversion and a 
corresponding inverse conversion. Thus the corresponding calculations can 
be performed faster than in the prior art solutions, such as systems based 
on FFT transformations. The method according to the present invention is 
characterized in what is set forth in the characterizing part of Claim 1 . The 
positioning system according to the present invention is characterized in 
what is set forth in the characterizing part of Claim 9. The receiver according 



to the first preferred embodiment of the present invention is characterized in 
what is set forth in the characterizing part of Claim 14. The electronic device 
according to the present invention is characterized in what is set forth in the 
characterizing part of Claim 22. The electronic device according to the 
second preferred embodiment of the present invention is further 
characterized in what is set forth in the characterizing part of Claim 29. 

The present invention provides considerable advantages as compared to 
the prior art methods and receivers. With the method according to the 
invention, the conversions can be performed faster and with lower power 
consumption as compared to, for example, FFT-based systems. Especially 
in location determination receivers, it is thus possible to make it faster for 
the receiver to become synchronized with the received signal. In addition, 
the real-time time-to-frequency and frequency-to-time conversion according 
to the invention can be implemented with a semiconductor circuit of a 
substantially smaller area than in the prior art solutions. The method 
according to the invention provides a substantial improvement of the 
synchronization speed as compared to the prior art methods. The total 
energy consumption of the receiver according to the invention can be kept 
reasonable, which makes the invention especially suitable for use in 
portable devices. Thus the location determination receiver can also be 
implemented in connection with a mobile station. 

In the following, the present invention will be described in more detail with 
reference to the accompanying drawings, in which 

Figure 1 is a simplified block diagram of a receiver, in which the method 
. according to the invention can be applied, 

Figure 2 shows the correlation step of a method according to a preferred 
embodiment of the invention as a simplified diagram, 

Figure 3 shows the analysjs step of a method according to a preferred 
embodiment of the invention as a simplified diagram, 



Figure 4 



shows the acquisition step of a method according to a preferred 
embodiment of the invention as a simplified diagram, 



Figure 5 shows an electronic device according to a preferred embodiment 
of the invention as a simplified block diagram, and 

Figure 6 illustrates the principle of a positioning system according to a 
preferred embodiment of the invention as a simplified diagram. 

In the following, the invention will be described using a GPS receiver 
according to Fig. 1 as an example. The properties of the Hartley transform 
and inverse Hartley transform applied in the method will be described at 
first. The Hartley transform can be presented with the formula: 

H ( f ) = X ) • [cos(27tfkt 0 / M) + sin(27rfkt 0 /M)]/M (1 ) 

and the inverse Hartley transform can be presented with the formula: 

*(0 = X "(^0 ) ■ [cos(2nkf 0 t I M) + sin(2nkf 0 t/ M)] (2) 
/f=i 

where 

x(k) are the measured samples in the time domain, 
M is the number of samples, 
to is the sampling interval, and 
f 0 is1/[(M-1)t 0 ] 

According to the formulas presented above, N 2 calculations are needed to 
calculate a discrete Hartley transform. However, this number can be 
reduced by means of permutations known from the theory of mathematical 
conversions and the so-called Butterfly process, whereby the conversions 
can be carried out with Nlog 2 N calculations. This conversion is then called a 
fast Hartley transform, FHT. On the basis of the properties of the Hartley 
transform, convolution in the time domain can be performed in the frequency 
domain by means of multiplications and additions as follows: 

*1 (0* MO = «i (f) ■ H 2e (f) + "i (-0 • H 2o (f) (3) 



where 
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x^t) and x 2 (t) are the signals on which convolution is performed, 
H,(f) and H 2 (f) are the x^t) and x 2 (t) Hartley transforms, 
respectively, and 

H 2e (f) and H 2o (f) are the even and odd transformation of the 
Hartley transform of the signal x 2 (t), respectively. 

If at least one signal of the convolution is odd or even, it follows that the 
convolution equation is simplified to the form H 1 (f)-H 2 (0- Tne above 
described convolution can be utilized when calculating, for example, the 
correlation between two signals. In order to calculate correlation by means 
of convolution, either of the signals is turned to its mirror image in the time 
domain, and the convolution is preferably calculated by using the above 
described Hartley transforms. This will be described later in this 
specification. 

The operation of a receiver 1 according to the first preferred embodiment of 
the invention will be discussed in the following. In the receiver 1, the 
received signal is preferably converted to an intermediate frequency in the 
converter block 2. At this stage, as is well known, the signal comprises two 
components, I and Q, between which there is a phase difference of 90°. 
These analog signal components, which are converted to the intermediate 
frequency, are digitized in the digitizing block 3 and conducted to the 
multiplier block 4. In the multiplier block 4, the I and Q components of the 
digitized signal are multiplied with the signal generated by the Numerically 
Controlled Oscillator (NCO) 5. This signal of the numerically controlled 
oscillator 5 is intended to correct the frequency deviation caused by the 
Doppler shift and the frequency error of the local oscillator (not shown) of 
the receiver 1. The signal formed by the multiplier block 4 is led to the 
synchronization block 6, in which the method according to the invention is 
applied. This synchronization block 6 thus tries to establish the code phase 
and frequency deviation of the satellite, which are to be used in functions 
after the synchronization. This will be discussed later in this specification. 
The scanning block 8 is controlled with the control block 7 during 
synchronization, and the scanning block 8 is used to regulate the frequency 
of the numerically controlled oscillator 5 when required. The control block 7 
controls the first switch 9 to switch either the signal formed by this scanning 
block 8 to the numerically controlled oscillator 5 during synchronization or 
the control signal formed by the code tracking block 1 1 to the numerically 



controlled oscillator 5 when synchronization has been achieved. A second 
switch 10 is used to control the operation of this tracking block 11. This 
tracking block 1 1 forms a part of a code phase locked loop and a carrier 
phase locked loop, which are known as such (not shown). 

After switching the operating voltages on or in a situation where the receiver 
1 has not been able to receive a signal of the GPS satellites for a long time, 
a two-dimensional search step is preferably carried out in the receiver 1 for 
each satellite to which signal reception is targeted. In this two-dimensional 
search step, the purpose is to determine the carrier frequency and code 
phase of each satellite. This carrier frequency is influenced by the Doppler 
shift caused by the movement of the satellite, the multipath propagation and 
the inaccuracies of the local oscillator of the receiver. The frequency 
inaccuracy can be rather high, as much as ± 6 kHz, in which case the 
receiver 1 must search from the frequency band of approx. 12 kHz in 
relation to the actual transmission frequency (L1 = 1575.42 MHz). The 
receiver 1 does not know the exact code phase, either, and hence the 
receiver must also determine the code phase from 1023 possible code 
phases. This results in a two-dimensional search process, in which one 
dimension is a frequency deviation in the band of 12 kHz, and the other 
dimension is a code phase out of 1023 different code phases. In a method 
according to a preferred embodiment of the invention it is possible to 
examine a frequency band of approx. 500 Hz at a time, and hence the 
method is repeated 24 times when required, to cover the whole frequency 
band of 12 kHz to be examined. Obviously the values used in this 
specification only serve as examples which clarify the invention, but do not 
constitute a limitation thereof. The invention can also be applied in other 
than GPS systems, in which case the frequency values, code phases and 
the number of codes can vary. 

The operation of a method according to the first preferred embodiment of 
the invention in the receiver 1 according to Fig. 1 will be described in the 
following. In order to start synchronization, the scanning block 8 sets the 
frequency of the numerically controlled oscillator 5 such that the receiver 1 
advantageously receives the lowest frequencies of the frequency band, 
1575.414 MHz— 1575.4145 MHz in this example. The receiver can also 
specify the starting frequency so that the receiver utilizes previously 
determined location information and/or almanac information, whereby the 
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location determination can be made even faster. Samples of the received 
signal are saved in the sample vector formation block 12 preferably as 
complex sample vectors p(1), p(2)...p(N), each of which includes 1023 
samples in this preferred embodiment. In this preferred embodiment, the 
sample saving frequency of the sample vector formation block 12 is 
essentially the same as the chipping rate of the chips, which is approx. 
1 023 000 samples per second. The sample vectors are continuous so that 
the next sample vector continues temporally after the previous sample 
vector, or the time interval between the last sample of the sample vector and 
the first sample of the next sample vector is essentially the same as the time 
interval between consecutive samples of the sample vector. These 1023 
samples thus correspond to a signal of 1 ms, and hence in a time-to- 
frequency transformation the frequency band is 1 kHz, which can be partly 
utilized. The sample vector formation step is denoted with the reference 
number 101 in Fig. 2 in this specification. 

The number of sample vectors is preferably N. In addition, the formation of 
sample vectors can be repeated. When the value of the number N of sample 
vectors p(1), p(2)...p(N) is determined in the GPS system, it must be taken 
into consideration that information has been modulated on to the signal at a 
bit rate of 50 bits/s as binary phase modulation. Another factor, which limits 
this number N of the sample vectors p(1), p(2)...p(N), is the frequency 
stability of the local oscillator of the receiver. 

In addition to the sample vector formation step, the synchronization method 
according to the invention also comprises a correlation step for forming a 
correlation function matrix. 

This correlation step can be partly performed already during sampling, or 
after N sample vectors p(1), p(2)...p(N) have been formed. If the correlation 
step is performed in such a manner, for instance, that a time-to-frequency 
transformation is calculated for each sample vector after it has been saved, 
preferably with a fast Hartley transform (FHT), the same time-to-frequency 
transformer can be used in all N sample vectors p(1), p(2)...p(N). But if the 
correlation step is performed after the saving of N sample vectors, a 
separate time-to-frequency transformer must be used for each sample 
vector, or the time-to-frequency transforms are performed for different 
sample vectors consecutively in the same time-to-frequency transformer. 
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Fig. 2 shows the correlation step of the method, in which a correlation 
function matrix C x .is formed of the sample vectors p(1), p(2)...p(N). 

A discrete Hartley transform 102, most suitably a fast Hartley transform, 
FHT, is performed on each sample vector p(1), p(2)...p(N). 

P(i) = FHT(p(/)) , where i=1 N (4) 

This is illustrated by the blocks FFT1, FFT2, .... FHTN in Fig. 2. The number 
of values used in calculations is preferably 1024, because then the discrete 
Hartley transform can be implemented in practical applications much more 
efficiently (with an FHT algorithm) than when 1023 values are used. One 
way of doing this is to add an extra null as the 1024th element. This has an 
insignificant effect on the transformation result. 

The reference codes r(x) corresponding to the C/A code of all the satellites 
of the GPS system are preferably saved in the receiver, whereby x refers to 
the satellite identifier and is, for example, in the range 1 — 30. It is not 
necessary to save the reference codes, but they can also be generated in 
the receiver. The reference code, which is selected or generated in the 
correlation step, is the reference code of the satellite transmitting the signal 
with which the receiver is being synchronized at any given time. The 
reference code is temporally inverted. This inverse reference code, which is 
denoted with r(x) in Fig. 2, is treated with a discrete Hartley transform 103, 

most suitably a fast Hartley transform: 

R(x) = FHT(f(x)) (5) 

The inverse reference code r(x) and/or its FHT transformation has often 

been saved in the memory media of the receiver in advance, but it can be 
formed of the reference code r(x) in connection with synchronization. 

A multiplication 104 between the Hartley transform result P(i) of each 
sample vector p(i) and the Hartley transform r(jc)of the inverse reference 
code R(x) is performed next in the correlation step. 
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M x (i)=R(x)-P(i) (6) 

Finally, an inverse Hartley transform 105 (IFHT) is performed on the results 
of these multiplications, resulting in a cross-correlation of the reference code 
5 r(x) and the received signal with all possible integer delays (1023). 



m x (i) = iFHT(M x (/)) 



(7) 
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This is based on the above mentioned fact that the Hartley transform of the 
convolution of the time domain signals corresponds to the multiplication of 
the Hartley transforms of the time domain signals. When an inverse 
reference code is also used, it is possible to perform a fast discrete-time 
correlation with the Hartley transformation. In this preferred example, the 
correlation result then comprises 1023 elements. The cross-correlation 
results m x (i) formed of these different sample vectors p(i) are used to form a 
correlation function matrix C x , in which the number of rows is the number N 
of sample vectors. 

It is obvious that instead of temporal inversion of the reference code, it is 
possible to form inverse sample vectors p(/') of the sample vectors p(i), in 

which case the reference code r(x) is used as direct and the sample vectors 
are used as inverse in the calculations presented above. In a preferred 
embodiment, it is not necessary to perform either of the inversions 
mentioned above, but the reference code r(x) and the sample vectors p(i) 
can be used directly. 

The rows of a correlation function matrix C x formed in the correlation step 
show the cross-correlation of the received signal and the reference code 
with different phase differences at the intervals of one millisecond. As a 
formula, the correlation function matrix can then be expressed as follows: 



(8) 



mx,kC0 
mx.k(2) 

m Xik (/v) 



At the next or analysis phase, the transpose 106 of the correlation function 
matrix C x is used; in the transpose the rows represent the signal samples in 

the time domain in the same manner as in a prior art correlator. Each row 
corresponds to a certain code phase difference between the received signal 
and the reference code. A Hartley transform 107 is performed on each row 
of the transpose of this correlation function matrix C x for forming a coherent 
search matrix A x , whereby a frequency analysis can be performed for 
determining the real frequency shift. This is illustrated by Fig. 3. 
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a x =fht(c t ) (9) 

In practical applications, it is not necessary to form a transposed matrix from 
the correlation function matrix in a separate step, but the elements of the 
saved correlation function matrix C x are read from the memory 16 (Fig. 5) in 
a different direction, preferably by columns. 

If necessary, an incoherent summing step can be performed in the 
synchronization block 6 to improve the signal-to-noise ratio. In order to 
implement the incoherent summing step, the sample vector formation step, 
the correlation step and the analysis step described above are repeated K 
times (not shown). This number of repetitions K is preferably selected such 
that the signal-to-noise ratio can be improved sufficiently, but in a 
reasonable time. 

The values of the elements a x (i,j) of the search matrix A x are examined in 
the acquisition step, and the objective is to find 110 a value, which exceeds 
the threshold value set in advance and is clearly higher than the other 
values. If a value like this is found, it expresses the code phase difference 
and frequency deviation, because it is probably a signal transmitted by a 
satellite. If the signal is not transmitted by a satellite, but consists of noise or 
other incidental spurious signals, no remarkable correlation peaks should 
arise. The code phase difference is seen from the row index of this highest 
value, and correspondingly the frequency deviation is seen from the column 
index. The example of Fig. 4 illustrates a situation where a value, which is 
clearly higher than the other values, has been found. If no such value is 
found in the search matrix A x , which means that a signal transmitted by the 
satellite searched for was probably not received in the frequency range 
examined, the frequency range to be examined is changed, and the above 
mentioned steps are carried out for forming a search matrix. With this 
method, the whole band of 6 kHz can be examined by repeating the steps 
described above as many times as required. 

When required, the steps described above can be repeated for the whole 
frequency range to be examined, and the search matrices formed at 
different times of repetition or only the possible peaks can be saved before 
searching for the highest correlation peak. Then the possibility of false 
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interpretations can be reduced in such a situation, for instance, that the 
threshold value has been set too small and a spurious signal can cause a 
false interpretation. 

5 After the right frequency deviation and code phase have been determined, 
the receiver can be set to the tracking mode. Reception of data is not 
successful with weaker signals, and data received via the mobile phone 
network, for example, must be used instead, which is a situation known as 
such. Range measurement is still possible with lower accuracy. The 
10 tracking mode is set by changing the switches (Fig. 1) to another position, 
whereby the received information is led to the tracking block 1 1 , in which a 
feedback is also formed for the fine adjustment of the frequency of the 
numerically controlled oscillator 5. 

15 For calculating the location, the receiver performs the signal reception 
preferably on the basis of signals received from at least four satellites. Then 
the synchronization described above is repeated when required for the 
signal of each satellite, and the code of the satellite with which 
synchronization is performed is selected as the reference sequence r(x) at 

20 each time. 

Most of the blocks needed in the implementation of the method according to 
the first preferred embodiment of the invention can be implemented in the 
digital signal processor (not shown), for example. Either hardware-based 
25 solutions or software implementations of the digital signal processing unit 
can be used in the performance of FHT transformations. In addition, a 
control element, preferably a microprocessor or the like can be used to 
control the operation of the receiver. 

30 Using a Hartley transform as the time-to-frequency and frequency-to-time 
transformation also provides the advantage that the same algorithm can be 
used in both transformations, which can be seen by comparing the formulas 
(1) and (2). In a hardware-based solution, the electronics required by the 
transformations can be implemented with either a circuit area reduced by 

35 about a half (an ASIC circuit, for instance) and half lower power 
consumption than a corresponding implementation with a Fourier transform 
and about the same processing time, or with a circuit area of the same size 
and the processing speed even doubled. 
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The transformation results H(f) of the Hartley transforms can also be used to 
calculate the amplitude spectrum and power spectrum of the signal. The 
amplitude spectrum can be calculated with the formulas: 

Xreal = H(f) + H i N ~ f) (real part) (1 0) 

Ximag = H{ f ) ~ H i N - 0 (imaginary part) (1 1 ) 

The power spectrum can be calculated with the formula: 

S = [H{ff + H(N-f) 2 y2 (12) 

An electronic device 24 according to a preferred embodiment of the 
invention, which comprises a mobile station and a location determination 
receiver, is shown in Fig. 5. The first antenna 13 is used for the reception of 
the signal transmitted by the positioning satellites. The received signal is led 
to the first radio part 14, in which the signal is down-converted and digitized. 
The first radio part comprises, among other things, the transformation block 
2, digitizing block 3 and multiplier block 4 of the receiver shown in Fig. 1 . 
The digitized signal, which at this stage preferably comprises the I and Q 
components, is led to the digital signal processing unit 15, in which, among 
other things, sample vectors are formed in the sample vector formation 
block 12. Samples are saved in the first memory media 16, which comprise, 
for instance, read/write memory and preferably also read memory and/or 
non-volatile read/write memory for saving the program code of the digital 
signal processing unit 15. In this embodiment, the signal processing unit 15 
also comprises operations of the synchronization block 6, for instance, such 
as the formation of the correlation function matrix C xk with, for example, 

Hartley transformers FFT1, FFT2, FHTN and/or matched filters. The 
digital signal processing unit 15 preferably also includes the formation of a 
coherent search matrix A x and the performance of an acquisition step. The 
digital signal processing unit 15 transmits information about the calculated 
phase difference and frequency deviation to the processor block 17, which 
comprises, for example, a microprocessor and I/O logic. The processor 
block 17 controls the scanning block 8 and the first switch 9. The tracking 
block 1 1 is preferably at least partly implemented as program instructions of 
the processor block. The second memory media 18 are used as the data 



memory and program memory of the processor block 17. It is obvious that 
the first memory media 16 and the second memory media 18 can also 
comprise common memory. Location determination information can be 
shown to the user on the display 1 9. 

Operations of the mobile station have also been implemented in the 
application software of the processor block 17. The display 19 can then be 
used for showing telephone call information in a manner known as such. 
With the keyboard 20, the user can control the location determination 
receiver and the mobile station. The coding and decoding of audio signals is 
performed by means of a codec 21. The radio part 22 and the second 
antenna 23 of the mobile station are also shown in Fig. 5. 

In the method according to the first preferred embodiment of the invention 
described above, location determination was performed in connection with 
the receiver. In a method according to another preferred embodiment of the 
invention, at least part of the location determination calculation is performed 
in connection with the data transfer network. A data transfer network can 
comprise a mobile communication network, for example. In that case, a 
calculation server 26 (Fig. 6) or a corresponding unit suitable for location 
determination calculation has been arranged in a data transfer connection in 
the data transfer network. One possible location for the calculation server is 
a base station 25 of the mobile communication network, but other alternative 
locations can also be applied in connection with the present invention. 

The operation of a method according to another preferred embodiment of 
the invention in an electronic device according to Fig. 5 and a system 
according to Fig. 6 will be described in the following. In this method, the 
operations performed are essentially the same as in the receiver 1 
according to the first preferred embodiment of the invention, and thus only 
the main features of the method will be described in the following. In order to 
start synchronization, the scanning block 8 sets the frequency of the 
numerically controlled oscillator 5 such that each time the frequencies 
wanted are received with the receiver 1. Samples of the signal being 
received are saved in the sample vector formation block 12, preferably as 
complex sample vectors p(1), p(2)...p(N). 
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Also in this embodiment, the correlation step can be partly performed 
already during sampling, or after N sample vectors p(1), p(2)...p(N) have 
been formed. If the correlation step is performed in such a manner, for 
instance, that a time-to-frequency transform is calculated for each sample 
vector after it has been saved, preferably with a fast Hartley transform 
(FHT), the same time-to-frequency transformer can be used in all N sample 
vectors p(1), p(2)...p(N). But if the correlation step is performed after the 
saving of N sample vectors, a separate time-to-frequency transformer must 
be used for each sample vector, or the time-to-frequency transforms are 
performed consecutively for different sample vectors in the same time-to- 
frequency transformer. In this embodiment, the correlation step is performed 
in the calculation server 26. Then the values of the samples of the sample 
vectors and the time data of the electronic device 24 are sent preferably by 
means of the radio part 22 of the mobile station to the mobile 
communication network. In the mobile communication network, these data 
are received in the base station 25 and transmitted to the calculation server 
26, where they are moved for processing to the block which performs the 
correlation step, depicted by the synchronization block 6 in Fig. 6. A 
correlation function matrix C x is formed of the sample vectors 

P0)> P(2)...p(N) in the manner specified in connection with the description of 
the method according to the first preferred embodiment of the invention. 

A discrete Hartley transform 102, most suitably a fast Hartley transform, is 
performed on each sample vector p(1), p(2)...p(N) in accordance with the 
formula (4) presented above in this specification. 

The reference codes r(x) of preferably all satellites of the GPS system (not 
shown), which correspond to the C/A code, have been saved in the 
calculation server 26, or they are generated in the calculation server 26. The 
reference code of the satellite to the signal of which the receiver is being 
synchronized is selected or generated each time in the correlation step. The 
reference code is temporally inverted. This inverse reference code, which is 
denoted by f(x) in Fig. 2, is treated with a discrete Hartley transform 103, 

most suitably a fast Hartley transform, in accordance with the formula (5) 
presented above in this specification. 
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A multiplication 104 between the Hartley transform result P(i) of each 
sample vector p(i) and the Hartley transform R(x) of the inverse reference 

code f(x) is performed next in the correlation step in accordance with the 
formula (6) presented above in this specification. 

Finally, an inverse Hartley transform 105 (IFHT) is performed on the results 
of these multiplications, resulting in a cross-correlation of the reference code 
r(x) and the received signal with all possible integer delays (1023) (formula 
7). 

In the analysis step, a Hartley transform is performed on each row of the 
transpose of this correlation function matrix C x for forming a coherent 
search matrix A x , whereby a frequency analysis can be performed for 
determining the real frequency shift. 

If necessary, an incoherent summing step can be performed in the 
synchronization block 6 to improve the signal-to-noise ratio. In order to 
implement the incoherent summing step, the sample vector formation step, 
the correlation step and the analysis step described above are repeated K 
times (not shown). This number of repetitions K is preferably selected such 
that the signal-to-noise ratio can be improved sufficiently, but in a 
reasonable time. 

The values of the elements A x of the search matrix a x (/,/) are examined in 
the acquisition step, and the objective is to find 110 a value, which is 
probably caused by a signal transmitted by a satellite. If such a value is 
found, the calculation server 26 transmits information of the code phase and 
the frequency shift to the electronic device 24. But if no such value is found 
in the search matrix A x , which means that a signal transmitted by the 
satellite searched for was probably not received in the frequency range 
examined, the frequency range to be examined is changed, and the above 
mentioned steps are carried out for forming a search matrix. In order to 
achieve this, the calculation server 26 transmits to the electronic device 24 
information that no signal of the satellite was detected in the frequency 
range being examined, whereupon the receiver 1 of the electronic device 
preferably changes the frequency range to be examined and starts 
performing the steps presented above in this changed frequency range. 
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With this method, the whole band of 6 kHz can be examined by repeating 
the steps described above as many times as required. 

When required, the steps described above can be repeated for the whole 
frequency range to be examined, and the search matrices formed in the 
calculation server 26 at different times of repetition or only the possible 
peaks can be saved before searching for the highest correlation peak. Then 
the possibility of false interpretations can be reduced in such a situation, for 
instance, that the threshold value has been set too small and a spurious 
signal can cause a false interpretation. 

The method according to another preferred embodiment of the invention is 
especially efficient because, among other reasons, there is a very small 
inaccuracy in the location readings of the calculation server 26 used for 
location determination, less than one millimeter in practice. This is due to 
the fact that the location determination server is located in a place known 
exactly (a base station, for example), and it usually has the Ephemeris data 
and almanac data of the satellites available each time. 

Thus one basic idea of network-based positioning is to implement all 
demanding calculations in a calculation server 26 or a corresponding device 
located in the data transfer network. The receiver 1 of an electronic device 
must then become locked to the signal of at least four satellites and add 
time information to their code signals for calculating pseudoranges. After 
this, the electronic device 24 transmits signals to the calculation server 26 
for calculating the actual location. 

In network-based positioning, the receiver 1 may consist, at the minimum, of 
only the equipment required for carrying out the chip and code phase 
measurement and for transmitting the measurement results via the network 
to the server. 

Although the invention was described above mainly in connection with the 
synchronization of the receiver, it is obvious that the invention can also be 
applied in connection with maintaining the synchronization. 
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The present invention is not limited to the above described embodiments 
only, but its details can be modified without departing from the scope 
defined by the attached claims. 
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Claims 

1. A method for performing the operations for synchronizing a receiver (1) 
with a received code-modulated spread spectrum signal, in which method at 
least one reference code (r(x)) is used, which reference code corresponds 
to a code used in the modulation, the frequency shift of the received signal 
and the code phase of the code used in the modulation are acquired, 
cfrwaclteirozed in that the method includes at least the following steps: 

- sample vector formation step, in which samples (101) are taken from the 
received signal for forming the sample vectors (p(1), p(2)...p(N)), 

- correlation _ step (102, 103, 104, 105), in which a first Hartley 
transform (R(x)) is formed on the basis of said reference code (r(x)), and 

a second Hartley transform (P(i)) is formed on the basis of each sample 
vector (p(1), p(2)...p(N)) ) a multiplication is performed between the first 
Hartley transform (R(x)) formed on the basis of said reference code and 

the second Hartley transform (P(i)) formed on the basis of each sample 
vector (p(1), p(2)...p(N)), an inverse Hartley transform is performed on 
each multiplication result (M x (i)), and 

- acquisition step (1 1 0), in which the frequency shift and code phase are 
acquired on the basis of the inverse Hartley transforms (m x (i)) of the 
multiplication results (M x (i)). 

2. A method according to Claim 1, clhiairacfleriized in that in the 

method the frequency range to be examined is specified, the specified 
frequency range is divided into two or more parts, whereupon said sample 
vector formation step and correlation step are performed on each part, and 
that the method also comprises an analysis step, in which values of the 
inverse Hartley transforms (m x (i)) of the multiplication results (M x (i)) are 
saved for forming a coherent search matrix (A x ), the acquisition step is 

performed after the examination of the specified frequency range, and the 
frequency shift and code phase are acquired on the basis of the highest 
value of the coherent search matrix ( A x ). 
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3. A method according to Claim 2, characterized in that a 
threshold value is specified in the method, and that quantity values of the 
elements (a x (/,/)) of the coherent search matrix (A x ) that exceed said 

threshold value are used in the acquisition step for acquiring the frequency 
shift and code phase. 

4. A method according to Claim 2 or 3, characterized in that in the 
method said sample vector formation step, correlation step and analysis 
step are repeated (108) for forming at least two coherent search matrixes 
(A x ), and that a summing step (109) is also performed in the method, in 
which summing step an incoherent search matrix (s x ) is formed by 
summing incoherently the values (a x (/,/)) of the equivalent elements of the 
coherent search matrix (A x ) formed at each time of repetition, and that said 
incoherent search matrix (s x ) is used in said acquisition step for acquiring 
the frequency shift and code phase. 

5- A method according to Claim 4, characterized in that in the 
method the frequency range to be examined is specified, the specified 
frequency range is divided into two or more parts, whereupon said sample 
vector formation step, correlation step, analysis step and summing step are 
performed on each part, and values of the elements of the incoherent 
search matrix (s x ) are saved, and that the acquisition step is performed 
after the examination of the specified frequency range, and the frequency 
shift and code phase are acquired on the basis of the highest value. 

6- A method according to Claim 4 or 5, characterized in that a 
threshold value is determined in the method, and that quantity values of the 
elements of the incoherent search matrix (s x ) that exceed said threshold 
value are used in the acquisition step for acquiring the frequency shift and 
code phase. 

7 - A method according to any one of the claims 1 to 6, 
characterized in that an inverse code(r(x)) corresponding to said 

reference code (r(x)) is used in the correlation step to form the first Hartley 
transform (R(x)). 
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8. A method according to any one of the claims 1 to 6, 
characterized in that an inverse code corresponding to each sample 
vector (p(1), p(2)...p(N)) is used in the correlation step to form the second 
Hartley transform ( P(/') ). 

9. A positioning system, which comprises at least a receiver (1), 
synchronization means (6) for performing synchronization operations to a 
received code-modulated spread spectrum signal, means (16) for using at 
least one reference code (r(x)) in connection with the synchronization, which 
reference code (r(x)) corresponds to a code used in the modulation, and 
means (15) for acquiring the frequency shift of the received signal and the 
code phase of the code used in the modulation, characterized in that the 
synchronization means comprise at least: 

- sample vector formation means (12, 16) for forming sample vectors 
(p(1), p(2)...p(N)) from the received signal, 

- correlation jjieans (6, FHT1, FHT2...FHTN) for forming a first Hartley 
transform (R(x)) on the basis of said reference code (r(x)), and for 

forming a second Hartley transform (P(i)) on the basis of each sample 
vector (p(1), p(2)...p(N)), means for performing a multiplication between 
the first Hartley transform (R(x)) formed on the basis of said reference 

code and the second Hartley transform (P(i)) formed on the basis of each 

sample vector (p(1), p(2)...p(N)), and means for performing an inverse 

Hartley transform on each multiplication result (M x (i)) for forming a 
correlation function matrix (C x ), 

- acquisition means (7, 17) for acquiring the frequency shift and code 
phase by using the values of the inverse Hartley transforms (m x (i)) of the 
multiplication results (M x (i)). 

10. A positioning system according to Claim 9, characterized in that 
the received signals are signals transmitted by satellites of the GPS system. 
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11. A positioning system according to Claim 9 or 10, which 
comprises a data transfer network, charsiefleirSzedl in that at least part of the 
synchronization means (6) are formed in connection with the data transfer 
network, and that a data transfer connection is arranged to be established 
between the data transfer network and the receiver (1 ). 

12. A positioning system according to Claim 11, chairaicteirnzedl in 
that the data transfer network comprises at least a mobile communication 
network. 

13. A positioning system according to Claim 9 or 10, chairaciteirizetil 
in that the synchronization means (6) are formed in the receiver (1). 

14. A receiver (1), which comprises at least synchronization means 
(6) for performing synchronization operations to a received code-modulated 
spread spectrum signal, said receiver (1) has means (16) for using at least 
one reference code (r(x)) in connection with the synchronization, the 
reference code (r(x)) corresponding to a code used in the modulation, and 
means (15) for acquiring the frequency shift of the received signal and the 
code phase of the code used in the modulation, cltoairaclteirozedl in that the 
synchronization means comprise at least: 

- sample vector formation means (12, 16) for forming sample vectors 
(p(1), p(2)...p(N)) from the received signal, 

- correlation jneans (6, FHT1, FHT2...FHTN) for forming a first Hartley 
transform (R(x)) on the basis of said reference code (r(x)), and for 

forming a second Hartley transform (P(i)) on the basis of each sample 
vector (p(1), p(2)...p(N)), means for performing a multiplication between 
the first Hartley transform (R(x)) formed on the basis of said reference 

code and the second Hartley transform (P(i)) formed on the basis of each 

sample vector (p(1), p(2)...p(N)), and means for performing an inverse 

Hartley transform on each multiplication result (M x (i)) for forming a 
correlation function matrix (C x ), 

- acquisition means (7, 17) for acquiring the frequency shift and code 
phase by using the values of the inverse Hartley transforms (m x (i)) of the 
multiplication results (M x (i)). 
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15. A receiver (1) according to Claim 14, characterized in that it 
also comprises at least 

- means (17) for specifying the frequency range to be examined, means (5, 
8, 9) for dividing the specified frequency range into two or more parts, 
whereupon the formation of the sample vectors and the formation of the 
correlation function matrix (C x ) are arranged to be performed for each 

part, 

- means for forming a coherent search matrix ( A x ), 

- means (16) for saving the values of the elements (a x (i,j)) of the coherent 
search matrix (A x ), and 

- means (15) for determining the frequency shift and code phase on the 
basis of the highest value of the coherent search matrix ( A x ). 

16. A receiver (1) according to Claim 14 or 15,. characterized in 
that it comprises means (1 5) for specifying the threshold value, and means 
(15) for comparing the values of said threshold value and the values of said 
coherent search matrix (A x ) for determining the frequency shift and code 

phase. 

17. A receiver (1) according to any one of the Claims 14, 15 or 16, 
characterized in that said formation of sample vectors, formation of a 
correlation function matrix (C x ) and formation of a coherent search matrix 
( A x ) are arranged to be repeated at least two times for forming a coherent 
search matrix (A xk ), and that the receiver (1) also comprises summing 
means for forming an incoherent search matrix (s x ) by summing the 
values (a xk (i,j)) of the equivalent elements of the coherent search 
matrix (A x ) formed at each time of repetition, and that said incoherent 
search matrix (s x ) is used in the acquisition step for determining the 
frequency shift and code phase. 

18. A receiver (1) according to Claim 14, characterized in that it 
also comprises at least 

- means (17) for specifying the frequency range to be examined, 

- means (5, 8, 9) for dividing the specified frequency range into two or more 
parts, whereupon the formation of the sample vectors and the formation 
of the correlation function matrix (C x ) are arranged to be performed for 
each part, 
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- means for forming a coherent search matrix ( A x ), 

- means for summing the coherent search matrix (A x ) to the incoherent 
search matrix (s x ), 

- means (16) for saving the values of the elements of the incoherent search 
matrix (s x ), and 

- means (15) for determining the frequency shift and code phase on the 
basis of the highest value. 

19. A receiver (1) according to Claim 18, characterized in that it 
comprises means (15) for specifying the threshold value, and means (15) for 
comparing the values of said threshold value and the values of said 
incoherent search matrix (s x ) for determining the frequency shift and code 

phase. 

20. A receiver (1) according to any one of the claims 14 to 19, 
characterized in that the correlation means comprise means (16) for 
forming a Hartley transform (R(x)) of the inverse code (r(x)) corresponding 

to said reference code (r(x)). 

21. A receiver (1) according to any one of the claims 14 to 19, 
characterized in that the correlation means comprise means (16) for 
forming a Hartley transform (P(/)) of the inverse code corresponding to 
each sample vector (p(1), p(2)...p(N)). 

22. An electronic device (24), which comprises at least a location 
determination receiver (1), synchronization means (6) for performing 
synchronization operations of the location determination receiver (1) to a 
transmitted code-modulated spread spectrum signal, and in which the 
location determination receiver (1) includes means (16) for using at least 
one reference code (r(x)) in connection with the synchronization, the 
reference code (r(x)) corresponds to a code used in the modulation, and the 
electronic device comprises means (15) for determining the frequency shift 
of the transmitted signal and the code phase of the code used in the 
modulation, characterized in that the synchronization means comprise at 
least: 

- sample vector formation means (12, 16) for forming sample vectors 
(p(1), p(2)...p(N)) from the received signal, 
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- correlation jneans (6, FHT1, FHT2...FHTN) for forming a first Hartley 
transform (R(x)) on the basis of said reference code (r(x)), and for 

forming a second Hartley transform (P(i)) on the basis of each sample 
vector (p(1), p(2)...p(N)), means for performing a multiplication between 
the first Hartley transform (R(x)) formed on the basis of said reference 

code and the second Hartley transform (P(i)) formed on the basis of each 
sample vector (p(1), p(2)...p(N)), and means for performing an inverse 
Hartley transform on each multiplication result (M x (i)) for forming a 
correlation function matrix (C x ), and 

- acquisition means (7, 17) for acquiring the frequency shift and code 
phase by using the values of the inverse Hartley transforms (m x (i)) of the 
multiplication results (M x (i)). 

23. An electronic device (24) according to Claim 22, which 
comprises means (14, 15, 17) for determining the location of the electronic 
device (24) and means (18) for saving the location information, 
cDnairaclteiriizedl in that the electronic device (24) also comprises means (17) 
for specifying the frequency range to be examined, and means (17, 22, 23) 
for selecting the starting frequency from said frequency range on the basis 
of the location information saved in the receiver (1). 

24. An electronic device (24) according to Claim 22 or 23, 
clharactorizedl in that the correlation means comprise means (16) for 
forming a Hartley transform (R(x)) of the inverse code (r(x)) corresponding 

to said reference code (r(x)). 

25. An electronic device (24) according to Claim 22 or 23, 
characterized in that the correlation means comprise means (16) for 
forming a Hartley transform (P(/)) of the inverse code corresponding to 
each sample vector (p(1), p(2)...p(N)). 

26. An electronic device (24) according to any one of the Claims 22 
to 25, characterized in that it also comprises means for performing data 
transfer operations. 

27. An electronic device (24) according to Claim 26, characterized 
in that it comprises means (22, 23) for establishing a data transfer 
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connection to a data transfer network, whereby the means for determining 
the location of the electronic device (24) comprise means (22, 23) for 
transmitting information needed in the location determination to the data 
transfer network, and means (22, 23) for retrieving information used in the 
location determination from the data transfer network, whereby at least part 
of the location determination operations are arranged to be performed 1n the 
data transfer network. 

28. An electronic device (24) according to Claim 26 or 27, 
cHnairactej-iized in that the means for performing data transfer operations 
comprise at least means for performing mobile station operations. 

29. An electronic device (24), which comprises at least a receiver (1) 
and which electronic device (24) is intended for use in connection with a 
positioning system, which comprises at least synchronization means (6) for 
performing synchronization operations of the receiver (1) to a transmitted 
code-modulated spread spectrum signal, means (16) for using at least one 
reference code (r(x)) in connection with the synchronization, the reference 
code (r(x)) corresponding to a code used in the modulation, means (15) for 
determining the frequency shift of the transmitted signal and the code phase 
of the code used in the modulation, and a data transfer network, 
cfhairacfteirozecfl in that the electronic device (24) also comprises at least 
sample vector formation means (12, 16) for forming sample 
vectors (p(1), p(2)...p(N)) from the received signal, and transmission means 
(22, 23) for transmitting the sample vectors (p(1), p(2)...p(N)) and time 
information to the data transfer network, and that the positioning system 
also comprises at least 

- correlation _means (6, FHT1, FHT2...FHTN) for forming a first Hartley 
transform (R(x)) on the basis of said reference code (r(x)), and for 

forming a second Hartley transform (P(i)) on the basis of each sample 
vector (p(1), p(2)...p(N)), means for performing a multiplication between 
the first Hartley transform (R(x)) formed on the basis of said reference 

code and the second Hartley transform (P(i)) formed on the basis of each 
sample vector (p(1), p(2)...p(N)), and means for performing an inverse 
Hartley transform on each multiplication result (M x (i)) for forming a 
correlation function matrix (C x ), and 
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- acquisition means (7, 17) for acquiring the frequency shift and code 
phase by using the values of the inverse Hartley transforms (m x (i)) of the 
multiplication results (M x (i)). 

30. An electronic device (24) according to Claim 29, chairactteirizedl 
in that it also comprises reception means (23, 22) for receiving information 
about the acquired frequency shift and code phase from the data transfer 
network. 

31 . An electronic device (24) according to Claim 29 or 30, which is 
intended for use in connection with a positioning system, and which also 
comprises at least means (26) for determining the location of an electronic 
device, clharacieirozed in that the electronic device (24) also comprises 
means (22, 23) for receiving location information from the data transfer 
network. 

32. An electronic device (24) according to Claim 29, 30 or 31, 
characterized in that the means for performing data transfer operations 
comprise at least means for performing mobile station operations. 



(57) Abstract 



The invention relates to a method for performing 
synchronization operations of a receiver (1) to a 
received code-modulated spread spectrum 
signal. At least one reference code (r(x)) is used 
in the method, which reference code 
corresponds to a code used in the modulation. 
The frequency shift of the received signal and 
the code phase of the code used in the 
modulation are determined in the method. The 
method also includes at least the following steps: 
sample vector formation step, in which samples 
(101) are taken from the received signal for 
forming the sample vectors (p(1), p(2)...p(N)), 
correlation step (102, 103, 104, 105), in which a 
first Hartley transform ( R(x) ) is formed on the 

basis of said reference code (r(x)), and a second 
Hartley transform (P(i)) is formed on the basis of 
each sample vector (p(1), p(2)...p(N)), a 
multiplication is performed between the first 
Hartley transform (R(x)) formed on the basis of 

said reference code and the second Hartley 
transform (P(i)) formed on the basis of each 
sample vector (p(1), p(2)...p(N)), an inverse 
Hartley transform is performed on each 
multiplication result (M x (i)), and 
acquisition step (110), in which the frequency 
shift and code phase are acquired on the basis 
of the inverse Hartley transforms (m x (i)) of the 
multiplication results (M x (i)). 
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